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PREFACE 


The  Civil  Engineering  Laboratory  (CEL) ,  under  sponsorship  of  the 
Naval  Facilities  Engineering  Command,  is  engaged  in  a  program  of  devel¬ 
oping  numerical  tools  for  analyzing  cable  systems.  As  part  of  this  work 
the  NRL-written  DESADE  program  has  been  used.  In  the  course  of  using 
DESADE  several  improvements  and  additions  have  been  made.  The  modifica¬ 
tions  are  for  the  most  part  minor  and  serve  to  increase  the  accuracy  of 
the  mathematical  modeling  and  to  add  user  conveniences. 

This  report  documents  these  modifications  in  the  form  of  a  revised 
users  manual.  (Much  of  the  text  of  the  original  manual  is  used  without 
modification.)  To  distinguish  the  updated  program  from  the  original  NRL 
version,  the  program  has  been  renamed  DECEL1  for  DESADE,  CEL  version  one 

CEL's  modifications  were  in  five  general  areas: 

1 .  Force  calculations 

2.  Current  field  definition 

3 .  Plotting  options 

4.  Miscellaneous  user  conveniences 

5.  Iteration  control 

Force  Calculations 


Drag  force  and  weight  of  cable  sections  covered  by  in-line  devices 
have  been  deleted.  (This  assumes  an  in-line  device  that  terminates  the 
cable  at  each  end  of  the  device;  floats  that  are  in-line  and  envelope  the 
cable  must  have  their  buoyancy  adjusted  by  the  weight  of  the  enveloped 
cable  segment.)  This  addition  has  a  relatively  minor  effect  except  when 
in-line  devices  envelope  a  significant  portion  of  a  cable. 

Calculation  of  tangential  drag  has  been  added  to  both  cables  and 
in-line  devices. 

Current  Field  Definition 


The  capability  to  accept  a  current  field  defined  by  2,  3  or  4 
current  meter  strings  of  up  to  25  meters  each  has  been  added.  Since 
much  current  meter  data  is  referenced  to  magnetic  north,  the  cable 
structure  can  be  referenced  to  the  current  field  by  specifying  the  angle 
between  the  structure's  X-axis  and  magnetic  north. 

The  specification  of  current  direction  has  been  changed  to  be 
consistent  with  oceanographic  terminology:  a  90°  current  flows  due 
east. 

An  option  has  been  added  to  specify  the  input  current  velocity 
units;  all  velocities  within  an  NDAT  case  must  be  the  same  units. 
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Plotting  Options 

Options  have  been  added  to  plot  the  current  field  defined  by  two  or 
more  current  meter  strings  and  to  plot  the  cable  structure  in  either  its 
deformed  or  undeformed  configuration.  Perspective  or  plan  or  elevation 
views  can  be  depicted. 

Miscellaneous  User  Conveniences 

Required  title  cards  have  been  added  for  the  cable  structure  source 
deck  and  for  each  parametric  case. 

The  specification  of  indexed  and  unindexed  devices  has  been  simpli¬ 
fied.  User  selected  devices  now  are  automatically  indexed,  in  order,  by 
their  location  along  a  cable. 

The  changes  allowed  in  a  parametric  study  deck  are  such  that  the 
physical  appearance  of  the  array  could  be  altered  in  the  parametric  case 
so  much  that  referencing  displacements  to  the  original  no-current  case 
is  illogical.  An  option  has  been  added  to  declare  any  parametric  case 
as  the  new  no-current  reference  case.  Displacements  are  printed  refer¬ 
enced  to  both  the  original  no-current  case  and  to  the  present  parametric 
reference  case.  The  reference  parametric  case  can  be  redefined  any 
number  of  times  since  the  original  no  current  case  is  retained  for  the 
duration  of  the  problem. 

The  error  detection  and  display  scheme  has  been  altered.  Errors 
that  are  detected  by  DECEL1  cause  the  entire  input  deck  to  be  listed. 
Then,  the  cards  with  errors  are  flagged  with  a  coded  error  number.  The 
coded  error  message  text  is  printed  below  the  input  card  listing.  All 
cards  are  scanned  for  errors;  however,  only  the  first  error  on  a  partic¬ 
ular  card  is  detected.  The  error  messages  are  identical  to  those  in  the 
original  DESADE  manual  except  that  only  the  portion  of  the  message 
applicable  to  the  card  type  is  printed. 

In  some  cases  it  is  desirable  to  be  able  to  punch  on  cards  the 
locations  of  particular  devices  for  input  to  other  programs.  An  option 
has  been  added  to  select,  based  on  device  weight,  the  devices  whose 
location  is  to  be  punched. 

Iteration  Control 

Under  some  circumstances  that  have  not  been  well  defined,  DECEL1 
may  either  fail  to  converge  or  converge  very  slowly.  To  protect  the 
user  from  high  execution  costs,  iteration  limits  have  been  imposed. 
Iterative  techniques  are  used  to  satisfy  the  imaginary  reaction  dis¬ 
placement  constraints  and  to  obtain  the  structure  shape.  Both  iteration 
processes  have  had  limits  imposed  because  both  have  been  the  cause  of 
excessively  high  execution  costs. 

These  modifications  have  added  to  the  capabilities  of  DECEL1  and 
have  made  it  a  more  useful  tool  for  the  cable  structure  analyst. 


INTRODUCTION 

DECEL1  is  a  Fortran  IV  program  for  computing  the  ocean  current- 
induced  static  deflections  of  undersea  structural  cable  arrays.  The 
solution  algorithm  is  the  Method  of  Imaginary  Reactions  (Refs  1,2) 


combined  with  the  method  of  successive  approximations  for  treating 
position  and  configuration  dependent  forces  (Refs  3,4). 

As  dimensioned,  the  program  can  handle  arbitrarily  configured 
arrays  of  up  to  22  cables.  The  cables  can  be  electromechanical,  wire 
rope,  or  synthetic.  Any  number  of  discrete  devices  (buoyancy  elements, 
current  meters,  tensiometers,  etc.)  may  be  incorporated  in  the  array. 

Certain  limitations  are  placed  on  the  structural  characteristics  of 
the  arrays  which  can  be  analyzed  by  this  program.  These  limitations  are 
as  follows: 

1.  No  cables  or  cable  segments  may  lie  on  the  ocean  floor. 

(No  surface  or  bottom  interaction  is  modeled,  thus  a 
cable  may  hang  below  the  lowest  anchor  point.) 

2.  The  dimensions  of  each  discrete  device  must  be  small  com* 
pared  to  overall  array  dimensions.  Thus,  for  example, 
the  application  of  this  program  to  the  analysis  of  an 
anchorage  for  a  submerged  submarine  is  not  valid. 

3.  All  parts  of  the  array  must  be  submerged.  Thus,  an  array 
containing  a  surface  buoy  cannot  be  validly  analyzed  using 
this  program.  (The  reason  for  this  is  that  a  surface  buoy 
generates  only  one  geometric  constraint  on  its  location.) 

An  exception  occurs  when  all  three  coordinates  of  a  device 
on  the  surface  can  be  specified  -  for  example,  the  coor¬ 
dinates  of  a  ship  handling  a  crown  line. 

Also,  certain  hydromechanical  assumptions  are  incorporated  in  the 
program  as  written.  These  are  as  follows: 

1 .  The  only  hydrodynamic  force  considered  to  be  acting  on 
the  discrete  devices  is  a  drag  force.  For  in-line  devices 
a  normal  drag  coefficient  and  a  tangential  drag  coeffi¬ 
cient  can  be  specified.  For  a  free  device,  only  the  nor¬ 
mal  drag  coefficient  is  to  be  specified.  Lift  forces  are 
neglected  as  being  small  compared  to  the  weight,  buoyancy, 
and  drag  forces  on  these  devices. 

2.  The  only  hydrodynamic  force  considered  to  be  acting  on  the 
cable  is  a  drag  force.  This  drag  force  consists  of  both 
normal  and  tangential  components  and  consequently  two 
drag  coefficients  (normal  and  tangential)  are  to  be  spe¬ 
cified  for  each  cable. 

3.  The  current  option  1  is  depth  dependent  in  magnitude,  uni¬ 
directional  and  horizontal.  The  current  option  2  permits 
direction  and  magnitude  variation  as  a  function  of  depth. 

The  current  option  3  also  permits  magnitude  and  direction 
variation  as  a  function  of  depth.  This  option  requires  two, 
three  or  four  "strings"  of  depth  dependent  velocity  data 
(where  the  strings  are  located  at  arbitrary  points).  These 
velocity  data  are  used  within  the  code  via  interpolation  or 
extrapolation  to  obtain  the  velocity  at  any  arbitrary  point 
in  space.  In  all  cases  the  vertical  component  of  velocity 
is  assumed  to  be  zero. 
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An  option  for  parametric  studies  (changes  in  weights,  diameters, 
cable  lengths,  anchor  locations,  etc.)  is  included  in  the  program. 

DECEL1  also  contains  a  series  of  error  checks  which  insure  that  all 
input  data  are  properly  formulated. 

Perspective  plotting  package  SSP  has  been  incorporated  into  DECELl 
to  plot  the  deformed  and/or  undeformed  configuration  of  the  array. 

The  program  uses  the  internally  generated  nodal  points  of  the 
cables  as  the  points  to  be  plotted. 

Scratch  files  3,  4  and  10  must  be  available  for  use  in  the  plotting 
package  SSP. 


Computer  Requirements 


As  written,  DECELl  should  be  compilable  on  most  Fortran  IV  compilers 
with  Boolean  algebra  capabilities.  Memory  requirements  for  the  program 
are  approximately  206,000  (octal)  words  in  single  precision.  Access  to 
one,  two,  or  three  magnetic  tape  units,  depending  on  the  I/O  options 
chosen,  is  required  by  the  program. 


Program  Operation 


The  overall  operation  of  the  program  DECELl  is  shown  in  the  flow 
diagram  in  Figure  1.  Numerical  examples  in  Appendix  A  and  the  source 
language  listing  of  the  program  is  given  in  Appendix  B. 


Array  Description 


Typical  cable  arrays  which  can  be  analyzed  using  DECELl  are  shown 
in  Figures  2a-2d.  The  figures  also  show  the  numbering  and  coordinate 
system  conventions  required  for  transmitting  the  array  geometry  to 
DECELl.  These  conventions  are  as  follows: 


1.  The  cables  comprising  the  array  must  be  numbered  consecu¬ 
tively  from  one  to  the  total  number  of  cables  in  the  array 
(Cl,  C2,  ...  in  Figure  2).  Each  cable  so  designated  must 
have  uniform  properties  (weight,  drag  coefficients,  diam¬ 
eter,  and  constitutive  relation)  along  its  length.  A 
change  in  property  also  requires  a  change  in  cable  number 
as  illustrated  in  Figure  2a. 

2.  The  termination  points  of  the  cables  in  the  array  are 
called  junctions.  A  junction  may  designate  an  anchor, 
the  intersection  point  of  two  or  more  cables,  or  the 
free  end  of  a  cable  such  as  illustrated  by  junction  J9 
in  Figure  2c.  The  junctions  must  also  be  numbered  con¬ 
secutively  from  one  to  the  total  number  of  junctions  in 
the  array  (Jl,  J2,  ...  in  Figure  2). 

3.  A  fixed,  right-handed  (X,  Y,  Z)  Cartesian  coordinate  sys¬ 
tem  must  be  chosen  to  describe  the  configuration  of  the 
array  in  space.  This  coordinate  system  is  called  the 
array-referenced  coordinate  system.  The  origin  of  the 
coordinate  system  can  be  arbitrarily  located.  The  Z 
axis  must  be  defined  parallel  to  the  direction  of  gravity 
and  increasing  upward.  All  distances  are  measured  in  feet. 
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4.  The  junction  numbers  corresponding  to  anchors  and  the  coor¬ 
dinate  of  the  anchors  must  be  tabulated  according  to  the 
scheme  illustrated  in  Table  1 . 


Table  1.  Anchor 

Tabulation 

Column  1 

Column  2 

Column  3 

Column  4 

Junction  number  of 

X  coordinate 

Y  coordinate 

Z  coordinate 

Anchor 

(ft) 

(ft) 

(ft) 

Anchors  are  defined  to  be  any  fixed  end  point  of  a  cable;  thus  an 
"anchor"  can  be  placed  on  the  bottom  or  at  the  surface  or  anywhere 
within  the  water  column. 

Reduction  to  a  Statically  Determinate  Array 

Before  an  arbitrary  cable  array  can  be  analyzed  by  DECEL1 ,  a  suffi¬ 
cient  number  of  cuts  must  be  made  in  the  array  to  reduce  it  to  a  stati¬ 
cally  determinate  structure.  The  effects  of  the  constraints  removed  by 
these  cuts  are  replaced  by  imaginary  and  equilibrating  reactions  (Refs 
1,2). 

The  number  of  cuts  required  to  reduce  a  cable  array  to  a  statically 
determinate  structure  is  determined  uniquely  from  the  relation 

number  of  cuts  =  number  of  cables  +  number  of  anchors 

-  number  of  junctions  (l) 


Certain  rules  must  be  adhered  to  as  the  required  cuts  are  made. 
These  rules  are  as  follows: 

1.  All  required  cuts  must  be  made  at  points  directly  adjacent 
to  array  junctions  -  that  is  at  end  points  of  the  cables 
comprising  the  array. 

2.  The  first  group  of  cuts  must  be  made  so  as  to  release  all 
but  one  cable  from  an  anchor. 

3.  The  remaining  cuts  (if  required)  are  made  within  the  array 
structure  and  must  be  located  so  as  not  to  break  the  array 
into  two  (or  more)  parts. 

4.  As  cuts  are  made,  each  new  cut  must  be  assigned  a  consecu¬ 
tive  junction  number,  continuing  from  the  last-used  junc¬ 
tion  number.  Also,  the  junction  number  (in  the  original 
array)  at  which  the  cut  is  made  must  be  tabulated. 

In  effect,  applying  rules  1-3  reduces  the  array  to  the  equivalent  of  a 
topological  tree.  As  the  name  implies,  this  is  a  continuous  structure 
containing  only  one  fixed  point  and  for  which  a  unique  (nonduplicative) 
path  exists  from  any  point  to  any  other  point. 


Examples  of  proper  reductions  to  statically  determinate  structures 
for  the  arrays  illustrated  in  Figures  2a-2d  are  shown  in  Figures  3a-3d, 
respectively.  In  each  of  these  figures,  the  left-hand  schematic  shows 
the  reduced  array  while  the  right-hand  schematic  depicts  the  topological 
tree  representation  of  the  reduced  array.  The  information  required  by 
rule  4,  which  represents  geometric  constraints  on  the  reduced  array,  is 
tabulated  below  the  left-hand  schematic  in  each  figure. 

Finally,  it  is  necessary  to  define  directions  of  increasing  arc 
length  along  the  cables  comprising  the  array.  The  tree  representation 
of  the  reduced  array  is  used  primarily  for  this  purpose.  These  direc¬ 
tions,  indicated  by  the  arrowheads  in  the  right-hand  schematic  of  each 
figure,  are  identified  by  starting  from  the  base  of  the  tree  and  climbing 
"up"  the  tree. 

Let  the  measure  of  arc  length  along  a  cable  be  denoted  by  s  which 
increases  from  zero  to  the  total  length  of  the  cable  L.  Then,  the 
required  information  on  increasing  directions  of  arc  length  can  be 
summarized  in  terms  of  array  junctions  as  shown  in  the  table  below  the 
right-hand  schematic  in  each  figure. 

Once  an  array  has  been  reduced  to  the  state  represented  by  Figure  3, 
it  is  amenable  to  analysis  by  the  program  DECEL1. 

Coordinate  Systems 

There  are  two  coordinate  systems  used  for  inputting  the  data  to 
DECEL1.  They  are  the  magnetically  aligned  reference  coordinate  system 
(N,  W,  Z)  and  the  array  or  laboratory  reference  coordinate  system  (X,  Y, 
Z) .  The  two  coordinate  systems  share  the  same  origin  and  the  same 
Z-axis;  Z  is  positive  upward.  Consequently,  the  two  coordinate  systems 
can  differ  from  one  another  by  an  arbitrary  angular  rotation  in  the 
horizontal  plane  which  is  denoted  by  $.  For  arbitrary  locations  on  the 
earth,  the  magnetic  axes  N-S  and  E-W  are  preestablished.  Consequently, 
the  angular  rotation  of  the  array-referenced  coordinate  system  is  refer¬ 
enced  to  the  magnetic  axes.  In  particular,  $  is  the  angle  between  the 
N-axis  and  the  positive  X-axis.  A  positive  rotation  of  the  X-axis  with 
respect  to  the  N-axis  is  in  the  clockwise  sense.  Figure  4  illustrates 
the  two  coordinate  systems . 

The  direction  of  positive  rotation  is  in  the  clockwise  sense  in  the 
magnetically  aligned  coordinate  system.  The  zero  degree  reference  is 
taken  as  the  N-axis. 

All  current  data  to  be  input  into  the  program  are  to  be  referenced 
with  respect  to  the  magnetic  coordinate  system.  For  example,  a  current 
having  an  inclination  of  0°  is  flowing  due  north;  a  current  with  270° 
inclination  is  flowing  due  west. 

When  current  field  input  option  3  is  selected,  it  is  necessary  to 
specify  the  locations  on  the  ocean  surface  (horizontal  plane)  of  the 
stations  where  current  data  have  been  gathered  (relative  to  the  N,W,Z 
coordinate  origin) .  Typical  examples  of  locations  of  such  stations 
might  be: 


100  ft  N  by  3000  ft  W  =  (  100,  3000) 


Cut  No. 


Junction  No. 
Assigned 
To  Cut 


Junction  No. 
At  Which 
Cut  Made 


1  6  1 

2  7  3 


(a) 


)6 


Junction  No. 
Assigned 

Cut  No.  To  Cut 


Junction  No. 
At  Which 
Cut  Made 


1 

2 

3 


7 

8 
9 


2 

3 

5 


(b) 


Cable  No. 

1 

2 

3 

4 


Junction  No. 
At  s  =  0 


4 
2 

5 
5 


Junction  No. 
At  s  =  L 


6 

5 

7 

4 


Junction  No. 
Cable  No.  At  s  *  0 


Junction  No. 
At  s  =  L 


1 

2 

3 

4 

5 

6 


1 

5 

6 
4 
6 
4 


4 

7 

8 

5 
9 

6 


Figure  3.  Examples  of  proper  reductions  to  statically  determinate  structures  for  the  arrays  shown  in 
Figures  2a  and  2b. 


CIO 


J8 


jn 


Junction  No. 
Assigned 

Cut  No.  To  Cut 


Junction  No. 
At  Which 
Cut  Made 


1  10  2 

2  11  3 

3  12  7 

4  13  7 


Junction  No. 
Cable  No.  At  s  =  0 


1  1 

2  5 

3  6 

4  4 

5  5 

6  6 

7  4 

8  8 

9  8 

10  7 


Junction  No. 
At  s  =  L 


10 

11 

12 

7 
13 

8 

5 

6 
9 


Junction  No. 
Assigned 

Cut  No.  To  Cut 


1 

2 

3 

4 

5 


8 

9 

10 

11 

12 


Junction  No. 
At  Which 
Cut  Made 


1 

2 

3 

4 

5 


Cable  No. 

1 
2 

3 

4 

5 

6 
7 

<d) 


Junction  No. 
At  s  =  0 

6 

6 

3 

7 

7 

7 

6 


Junction  No. 
At  s  =  L 

8 

9 

6 

10 

11 

12 

7 


Figure  3.  Examples  of  proper  reductions  to  statically  determinate  structures  for  the  arrays  shown  in 
Figures  2c  and  2d. 
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In  the  magnetic  coordinate  system  the  N,W  axes  are  positive.  (N,W)  form  a  coordinate 
pair  of  a  right  handed  coordinate  system.  The  positive  direction  of  rotation  0  in  the 
(N,W)  coordinate  system  is  clockwise  with  the  zero  degree  reference  being  the  N  -  axis. 
The  array  coordinate  system  (X,Y)  system.  The  relationship  between  the  two  systems 
is  given  above. 


Figure  4.  Description  and  relationship  between  the  magnetic  and  array  referenced  coordinate  systems. 


Above  (on  the  right),  is  an  example  of  how  such  locations  should  be 
interpreted  for  input  preparation  as  coordinate  pairs.  Note  that  the 
north-west  quadrant  is  the  first  quadrant  of  a  right-handed  coordinate 
system.  Consequently,  north  and  west  are  positive  axes  while  south  and 
east  are  the  negative  axes.  When  current  field  input  option  3  is  selected, 
the  current  data  are  printed  out  in  both  the  magnetic  and  array- referenced 
coordinate  systems. 

The  specification  of  the  anchor  locations  is  accomplished  with 
respect  to  the  array-referenced  (X,  Y,  Z)  coordinate  system.  All  internal 
calculations  within  DECEL1  are  performed  with  respect  to  the  array- 
referenced  coordinate  system.  The  (X,  Y,  Z)  coordinate  system  is  also 

right-handed.  In  the  horizontal  plane  in  this  frame  of  reference  the 

direction  of  positive  rotation  is  counter-clockwise  with  the  positive 
X-axis  being  the  zero  degree  reference. 

Directions  of  Positive  Rotation 

Reference  System  0°  Reference  Positive  Rotation 

Magnetic  N-axis  clockwise 

Array  X-axis  counter-clockwise 


Current  Field  -  Input  Option  0 

There  is  an  input  option  for  determining  the  static  deflections  of 
a  cable  system  immersed  in  a  currentless  environment.  This  is  current 
field  input  option  0.  This  option  is  automatically  processed  when  any 
of  the  subsequently  described  current  field  input  options  are  exercised. 

It  can  be  exercised  independently. 

Current  Field  -  Input  Option  1 

The  current  field  option  1  is  taken  to  be  unidirectional  and  hori¬ 
zontal,  though  depth  dependent  in  magnitude.  The  direction  of  the 
current  is  specified  from  the  magnetic  north  axis  by  using  the  ANG  card. 
[Within  the  code,  the  current  is  referenced  to  the  array  coordinate 
system  for  calculation  purposes.  Thus,  if  the  direction  of  the  flow 
with  respect  to  the  X-axis  is  denoted  by  6,  the  current  field  is  expressed 
by: 


Vj  =  V(Z)(ej  cos  0  +  e2  sin  0) 


Here,  e.,  e.  are  unit  base  vectors  with  respect  to  the  X-  and  Y-  axes, 
respectively.]  V(Z)  specifies  the  magnitude  of  the  current  as  a  function 
of  depth,  Z.  This  functional  relationship  must  be  tabulated  as  in 
Table  2. 


12 


Table  2.  Current  Tabulation  -  Option  1 


Column  1  Column  2 

Z  coordinate  (ft)  V(Z)  knots 


Up  to  25  rows  are  permitted  in  Table  2.  At  least  one  of  the  Z-coordinates 
in  Table  2  must  be  less  than  or  equal  to  the  Z-coordinate  of  the  lowest 
anchor.  A  sorting  scheme  is  invoked  within  the  program  that  arranges 
Table  2  data  according  to  ascending  values  of  Z.  Such  sorting  is  neces¬ 
sary  within  the  program  for  the  determination  of  the  current  of  an 
arbitrary  depth.  This  is  accomplished  via  a  linear  interpolation  of  the 
currents  using  the  data  at  the  two  depths  encompassing  the  depth  of 
interest. 

Figure  5  illustrates  the  magnetic-referenced  coordinate  system  and 
the  linear  interpolation  of  the  velocity  between  given  data  points. 

For  current  field  input  option  1,  the  directionality  of  the  current 
is  a  constant  for  the  entire  velocity  field.  Also,  for  any  point  (X,  Y) 
on  the  plane  (Z  =  constant),  the  value  of  velocity  is  invariant. 

Current  Field  -  Input  Option  2 

Current  field  input  option  2  is  a  slight  generalization  of  the 
current  field  input  option  1.  The  generalization  involves  allowing  the 
current  direction  8  to  vary  with  depth.  The  input  specification  for 
this  option  requires  depth  Z,  current  magnitude  and  current  direction 
from  the  magnetic  north  axis.  The  positive  direction  of  rotation  is 
clockwise  in  the  magnetically-aligned  reference  coordinate  system  and  a 
zero  degree  current  flows  due  north.  The  current  specification  for 
option  2  requires  a  tabulation  as  in  Table  3. 


Table  3.  Current  Tabulation  -  Option  2 
Column  1  Column  2  Column  3 


Z-coordinate  V(Z)  magnitude,  knots  8(Z)  direction,  degrees 

from  magnetic  N-axis, 
positive  in  clockwise 
sense 


Up  to  25  rows  are  permitted  in  Table  3.  The  velocity  sorting  scheme  is 
invoked  by  the  program.  One  entry  for  this  data  must  correspond  to  a 
depth  less  than  or  equal  to  the  Z-coordinate  of  the  lowest  anchor.  The 
above  current  magnitude  and  direction  data  are  used  to  generate  velocity 
components  along  the  (X,  Y)  axes.  To  obtain  the  velocity  at  an  arbitrary 
depth,  linear  interpolation  is  performed  using  the  (X,  Y)  velocity 
components  above  and  below  the  depth  of  interest;  the  angular  direction 
of  the  current  is  also  found  by  interpolation. 
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For  current  option  2,  the  value  of  velocity  at  any  (X,  Y)  point  on 
the  plane  Z  =  constant  is  invariant. 

An  ANG  card  can  be  used  to  rotate  the  entire  current  profile  in  the 
sane  manner  that  a  unidirectional  current  profile  is  rotated. 


Current  Field  -  Input  Option  3 


The  current  field  input  option  3  employs  a  linear  interpolation/ 
extrapolation  scheme  on  current  (magnitude  and  direction)  data  to  give 
velocity  variation  in  horizontal  planes.  When  this  option  is  invoked 
the  program  expects  current  data  from  2,  3  or  4  current  meter  strings. 

Up  to  25  measurements  of  current  data  can  be  contained  on  each  string. 

For  each  string  the  first  data  entry  must  be  at  a  depth  equal  to  or  less 
than  the  lowest  anchor.  The  depth  dependent  data  on  the  various  strings 
do  not  have  to  correspond  one  to  the  other.  That  is,  if  on  one  string 
current  data  are  obtained  at  Zj,  Z^,  Z„,  ...  Z  ,  then  on  the  other 
strings  data  can  be  collected  at  completely  di?ferent  depths.  This 
current  field  input  option  has  been  designed  to  treat  mildly  varying 
current  fields.  It  does  not  provide  acceptable  results  for  eddy  currents 
or  reverse  shear  flows.  The  interpolation/extrapolation  scheme  to 
obtain  the  velocity  at  an  arbitrary  point  in  the  field  works  as  follows: 


Case  of  2  Strings.  The  field  point  (X  ,  Y  ,  Z  )  is  determined 
where  the  velocity  is  to  be  found.  On  the  plane  Z=2  there  are  three 
points  of  interest.  They  are  (Xq,  Yq,  Zq),  (X^,  Y^Z  )  and  (X„,  Y2, 

Z  )  where  the  last  two  points  denotethe  location  on  tBe  plane  Z=Z  at 
tBe  meter  strings.  At  each  of  these  locations  current  magnitude  and 
direction  data  are  determined  by  a  straight  line  interpolation  of  the 
data  values  given  directly  above  and  below  the  Z  elevation.  Figure  6 
illustrates  the  situation  after  this  linear  interpolation.  There,  (V  , 
0j)  and  (V2>  02)  are  obtained  from  the  vertical  interpolation  just 
mentioned. 

A  line  normal  to  the  line  connecting  (Xj,  Y^)  and  (X2,  Y2)  through 
the  field  point  (X  ,  Y  )  is  constructed  to  locate  the  interpolation 
point  (X.,,  Y^).  TBe  pSint  (X~,  Y-)  lies  on  the  line  joining  (X.,  Y.) 
and  (X„,  Y2)T  Linear  interpolations  are  now  performed  to  find  v_  and  0. 
from  tne  corresponding  values  at  (Xj,  Yj)  and  (X2>  Y_).  The  current 
(V  ,  03)  are  assumed  to  prevail  at  all  points  along  the  normal  line. 
Hence,  (V_,  0.)  provide  the  description  of  the  current  at  the  field 
point  (XQf  Yj. 


Case  of  3  Strings.  Suppose  the  field  point  is  (X  ,  Y  ,  Z  ).  The 
current  magnitude  and  direction  are  determined  by  linear  vertical  inter¬ 
polation  at  the  three  string  locations  (X..,  Y  ,  Z  ),  (X»,  Y_,  Z  )  and 
(X_,  Y~,  Z  ).  The  corresponding  velocity1data  obtained  at  these  locations 
are  (Vj,  0°),  (V2,  0.)  and  (V_,  0_).  At  the  string  location  on  the 
horizontal  planez(Z=z  )  V.,  V2  ana  V_  are  amplitudes  through  which  a 
velocity  plane  can  be  passed.  (A  unique  plane  can  be  passed  through  any 
three  non-colinear  points.)  The  velocity  V  at  the  field  point  (X  ,  Y  , 

Z  )  then  can  be  obtained  by  determining  theamplitude  on  the  velocity  ° 
p?ane  corresponding  to  the  position  (X  ,  Y  ,  Z  ) .  The  procedure  is 
exactly  the  same  for  the  determination°of  Bhe  Current  direction  6  .  For 
this  quantity  a  direction  plane  has  to  be  established  using  the  direction 
amplitudes  0j,  0^  and  0^. 
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Line  0  -  3  is  orthogonal  to  line  1  -  2.  The  velocities  (Vj,  9 j)  and  (V2,  9^)  are  obtained  by 
linearly  interpolating  string  1  and  string  2  data,  respectively.  The  velocity  of  any  point  along 
line  1  -  2  can  be  obtained  by  linear  interpolation/extrapolation  of  the  velocities  (Vj,  9 j), 

(V2,  82)-  The  assumption  is  made  that  the  velocity  (Vj,  9 j)  of  an  arbitrary  point  (Xj,  Yj) 
on  line  1  -  2  is  propagated  along  the  orthogonal  trajectory  to  that  line  at  that  point. 


Figure  6.  Interpolation/extrapolation  in  the  case  of  2-string  data. 
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Case  of  4  Strings.  The  field  point  is  (X  ,  Y  ,  Z  ) .  The  string 
points  are  (X  ,  Y^  Z  ),  (X2,  Y2>  Zq)  ,  (X3>  Y3?  zj  an8  (X^  Y  Zq). 
Through  any  three  of  the  string  points  the  procedure  used  in  tne  preced¬ 
ing  case  can  be  repeated  exactly.  There  are  four  independent  ways  that 
three  points  can  be  selected  from  four  points,  namely:  (1,  2,  3);  (1, 

3,  4);  (1,  2,  4);  (2,  3,  4).  The  current  is  determined  at  the  field 
point  by  averaging  the  results  obtained  by  applying  the  case  of  three 
strings  four  times  to  the  four  combinations  of  points  just  indicated. 


TANGENTIAL  DRAG 

The  tangential  drag  computation  on  the  cable  is  straightforward  and 
is  based  on  the  expression 

"TD0  *  5  0  >"»<•«>'  CDTVTMVl(’j) 


/X\ 

where  AF  I  Y  I  =  tangential  drag  force  component 

TD  \  Z  /  per  unit  length  in  the  direction  X,  Y  or  Z 

/ X\ 

V  I  Y  1  =  tangential  velocity  component  along  X,  Y  or  Z 

T  \  Z  /  axis 

VTM  =  \^T(X)  +  VT(Y)  *  VT(Z)*  ma«nitude  of 

tangential  velocity 

Cjyj.  =  tangential  drag  coefficient 

p  =  density  of  fluid 

D  =  cable  diameter 

e  =  strain  developed  by  cable  tension 

The  above  calculation  forms  a  part  of  the  hydrodynamic  cable  drag  force. 
The  hydrodynamic  force  is  comprised  of  a  normal  drag  and  a  tangential 
drag. 

CABLE  DEVICES 

In  FUNCTION  EFORCE  (I,  M,  N)  the  tangential  drag  contribution  to 
in-line  devices  is  included  in  the  drag  calculation.  It  is  required 
therefore  to  supply  a  tangential  drag  coefficient  for  each  in-line 
device  as  input. 

For  in-line  cable  devices  the  weight  of  the  cable  and  the  hydrody¬ 
namic  drag  over  a  portion  of  the  cable  covered  by  an  in-line  device  is 
deleted.  Cable  devices  are  referred  to  as  in-line  or  free  devices  (two 
types).  All  cable  devices  are  numbered  automatically  within  the  program. 
The  numbering  is  accomplished  sequentially  with  respect  to  cable  number 


W.1! 


■59* 


and  location  of  the  device  away  from  its  s=o  end.  That  is,  if  on  cable 
1  there  are  three  cable  devices  located  at  s=s  ,  s=s2>  and  s=Sg  with  a 
s;j>s2>s1,  then  cable  device  indices  1,  2,  3  refer  to  aevices  located  on 
cable  1  at  s=Si,s2,Sg,  respectively. 

REFERENCE  CONFIGURATION 

An  option  exists  to  specify  a  particular  parametric  study  case  as  a 
new  no  current  reference  configuration.  Displacements  with  respect  to 
the  new  no  current  reference  configuration  as  well  as  the  displacement 
from  the  original  no  current  reference  configuration  are  part  of  the 
output. 


ITERATION  CONTROL  AND  EXECUTION  ERROR  MESSAGES 

To  protect  the  user  from  excessively  high  execution  costs  coupled 
with  the  risk  of  not  receiving  any  output,  iteration  limits  have  been 
added  to  DECEL1.  Two  separate  iteration  processes  are  involved:  one 
attempts  to  satisfy  the  displacement  constraints,  where  the  cuts  are 
made,  imposed  by  the  COMP  card;  where  the  other  determines  the  structure’s 
shape.  Both  iteration  processes  have  been  observed  on  rare  occasions  to 
converge  very  slowly,  if  at  all. 

The  first  iteration  process  deals  with  the  imaginary  reactions. 

The  associated  displacement  errors  have  been  observed  to  be  large  and  to 
change  slowly  in  some  slowly  convergent  problems.  An  arbitrary  definition 
of  slow  convergence  coupled  with  a  maximum  iteration  limit  is  used  to 
terminate  execution  of  a  particular  parametric  case.  Slow  convergence 
has  been  defined  to  occur  when,  after  half  of  the  iteration  limit  is 
achieved,  the  displacement  error  is  large  (100  times  the  COMP  value)  and 
the  error  is  changing  by  less  than  the  COMP  value  per  iteration.  Itera¬ 
tion  is  terminated  due  to  slow  convergence  with  the  message:  SLOW 
CONVERGENCE  AFTER  XXX  ITERATIONS.  The  previous  and  present  displacement 
error  values  are  printed  to  aid  in  determining  how  closely  the  solution 
has  converged.  If  only  one  of  the  slow  convergence  criteria  is  satisfied 
iteration  will  continue  until  convergence  occurs  or  the  iteration  limit 
is  reached.  Iteration  is  terminated  due  to  reaching  the  iteration  limit 
with  the  message:  NO  CONVERGENCE  IN  XXX  ITERATIONS  and  the  present  and 
previous  displacement  error  values  are  printed. 

The  second  iteration  process  to  determine  the  structure's  shape 
occurs  after  the  imaginary  reaction  iteration  has  been  successful.  This 
process  has  been  observed  on  some  occasions  to  oscillate  about  the 
correct  solution.  The  normal  iteration  process  is  allowed  to  continue 
until  half  of  the  nodal  position  iteration  limit  is  reached.  At  this 
point,  the  solution  is  arbitrarily  assumed  to  be  oscillating  and  a 
half-step  iteration  scheme  is  imposed  (the  iteration  process  is  unchanged 
except  that  each  node  is  allowed  to  move  only  half  as  far  as  calculated). 
This  half-step  technique  continues  until  the  iteration  limit  is  reached 
or  convergence  occurs.  Execution  of  the  parametric  case  is  terminated 
with  the  message:  PROGRAM  DID  NOT  CONVERGE  AFTER  XXX  ITERATIONS,  PARA¬ 
METRIC  CASE  TERMINATED.  The  displacement  error  is  printed  and  the  full 
output  for  the  case  is  printed  preceded  by  the  message:  APPROXIMATE 
RESULTS  PRINTED.  The  user  must  judge  if  an  adequate  displacement  error 
has  been  achieved  for  the  results  to  be  meaningful. 
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What  To  Do  If  .  .  . 


When  execution  has  been  terminated  by  reaching  the  imaginary  reaction 
iteration  limit,  several  options  are  still  available.  The  iteration 
limit  (field  4  of  the  COMP  card)  can  be  increased;  however,  this  probably 
will  have  to  be  done  in  conjunction  with  an  increase  in  the  COMP  value. 

The  COMP  value  should  still  be  kept  in  a  range  that  will  be  judged  to 
produce  usable  results.  The  error  values  printed  upon  execution  termina¬ 
tion  will  give  an  idea  of  how  much  the  COMP  value  will  have  to  be  increased. 

In  some  cases,  the  desired  current  is  too  strong  to  apply  in  one 
step;  strong  currents  have  caused  the  imaginary  reaction  iteration  limit 
to  be  reached.  To  handle  this  and  other  sensitive  cases,  an  option  has 
been  added  to  apply  the  current  in  increments  (field  11  of  the  NDAT 
card).  Using  this  option,  the  program  iterates  to  a  solution  for  the 
first  current  increment  as  if  this  were  the  total  current  to  be  applied. 

Then  the  current  is  incremented  and  the  process  is  repeated.  An  option 
exists  to  print  the  solution  for  each  increment  of  the  current  (see 
field  11  of  the  NDAT  card).  This  technique  has  been  highly  effective  in 
obtaining  solutions  in  high  currents  up  to  10  knots. 


DECEL1  USER  EXPERIENCE 

In  the  course  of  using  DECEL1  a  wide  variety  of  structures  in 
various  current  profiles  have  been  analyzed.  Some  results  have  been 
compared  with  experimental  data  with  close  agreement  (Ref  7).  Other 
cases  have  been  cross-checked  against  other  computer  models  again  with 
good  agreement.  The  program  is  easy  to  use  to  the  point  that  a  user 
with  a  general  engineering  background  can  use  the  manual  to  formulate 
the  required  input  and  expect  to  receive  back  usable  results  once  obvious 
input  errors  are  corrected.  The  ease  of  use  is  an  exceptional  attribute 
of  the  program. 

In  the  course  of  making  modifications  to  the  program  a  great  deal 
of  insight  has  been  gained  regarding  the  mathematical  modeling.  Many  of 
the  details  are  of  no  importance  to  the  user.  However  there  are  several 
characteristics  that  the  user  should  be  aware  of;  these  deal  with  the 
internal  distribution  of  loads,  non-convergent  problems,  and  the  conver¬ 
gence  parameter. 

Load  Distribution 


Loads  (weight  and  drag)  are  generated  from  three  sources:  cables, 
cable  devices,  and  junction  devices.  Cable  forces  due  to  drag  and 
weight  are  calculated  for  the  cable  length  between  two  nodes  [length  = 
(total  cable  length)/(no.  of  elements  in  the  cable)].  (See  CAB  card.) 
The  forces  are  assigned  to  the  M-l  node  as  shown  below: 


load 


Clearly,  the  smaller  the  distance  between  the  nodes,  the  more  accurate 
the  mathematical  representation. 

Loads  due  to  devices  on  the  cable  (DCAB  card)  are  modeled  in  a 
similar  way  except  that  the  device  is  physically  located  at  a  specific 
point  along  the  cable  between  the  two  nodes.  The  DCAB  loads  are  modeled 
as  shown  below: 


load 

t - 


•- - 

- - . 

M-l 

i 

DCAB  device 


This  is  a  good  representation  when  the  device  is  actually  near  node 
(M-l).  However,  when  the  device  is  near  node  M  and  the  distance  between 
the  nodes  is  large,  this  may  be  a  poor  representation.  This  is  not  to 
say  the  computer  results  will  be  grossly  in  error;  this  depends  upon  the 
size  and  number  of  DCAB  devices.  If  the  devices  are  few  and/or  small 
such  that  the  system  is  cable-load  dominated,  the  results  should  be 
acceptable. 

Probably  the  most  serious  error  will  arise  when  there  are  few 
elements  and  the  DCAB  devices  produce  large  static  or  drag  loads.  In 
this  case  large  loads  may  be  assigned  to  locations  quite  distant  from 
their  actual  point  of  application. 

Loads  due  to  DJNC  devices  are  most  accurately  modeled.  The  loads 
are  assigned  exactly  where  the  device  is  located. 

It  may  appear  that  a  better  modeling  of  cable  and  DCAB  loads  could 
be  achieved  by  proportioning  the  loads  between  the  appropriate  adjacent 
nodes.  However,  this  single  change  is  not  consistent  with  other  internal 
force  accounting  schemes.  To  implement  force  proportioning,  a  major 
re-write  of  the  program  would  be  required;  this  has  not  been  done. 

As  the  program  stands,  the  majority  of  problems  can  be  analyzed 
with  wholly  satisfactory  results  even  ignoring  the  way  loads  are  distri¬ 
buted.  Where  the  load  distribution  scheme  is  judged  by  the  user  to  pose 
a  problem  in  obtaining  a  satisfactory  model,  the  following  suggestions 
are  presented. 

1.  Model  all  cables  with  as  many  elements  as  possible  (50).  This 
will  make  the  distance  between  nodes  least  so  that  DCAB  loads  will  be 
assigned  as  close  to  the  actual  device  location  as  possible. 

2.  Model  large  bodies  as  DJNC  devices  -  loads  will  be  assigned 
exactly  where  the  device  is  located. 

3.  Model  long  cables  as  a  series  of  shorter  series-connected 
cables.  This  will  aid  in  making  elements  short. 
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Nog-Convergent  Problems 

DECEL1  treats  statically  determinate  and  statically  indeterminate 
problems  differently.  The  indeterminate  problems  are  solved  using  the 
imaginary  reaction  technique  (Ref  1).  This  is  a  powerful  method  and 
usually  converges  rapidly  to  a  solution.  (Reference  1  shows  that  the 
method  is  unconditionally  convergent.)  However,  at  least  one  case  has 
been  encountered  where  the  method  is  at  best  very  slowly  convergent. 
The  case  involves  a  nearly  neutrally  buoyant  cable  hung  in  a  catenary 
between  two  fixed  points  with  a  strong  current  acting  perpendicular  to 
the  plane  of  the  catenary: 


+Z 


current 

profile 


The  positions  attained  by  the  cable  during  the  iteration  procedure  are 
sketched  below. 


1 


The  cable  is  initially  driven  from  the  no-current  position  (1)  to  a 
nearly  horizontal  position  (2).  Here  there  is  low-  or  no-normal  drag  to 
support  the  cable  so  in  the  next  iteration  it  falls  to  position  (3) 
(which  may  be  coincident  with  position  (1)).  Internally,  the  iterations 
appear  to  oscillate  between  positions  (2)  and  (3),  and  the  solution 
fails  to  converge.  Convergence  can  be  achieved  by  reducing  the  current 
or  by  increasing  the  weight  of  the  cable. 

Statically  determinate  cases  bypass  the  imaginary  reaction  routine, 
thus  convergence  is  not  as  certain.  There  are  combinations  of  current 
and  geometry  that  interact  in  such  a  way  as  to  appear  to  be  unstable. 

An  example  it.  a  lightweight  anchor  suspended  from  the  surface  in  a 
current  that  increases  with  depth. 
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From  the  no-current  position  (1),  the  first  iteration  is  to  position  (2) 
where  the  normal  drag  force  is  reduced  both  because  of  the  orientation 
change  (toward  horizontal)  and  the  reduced  current  in  the  new  position. 
The  reduced  loading  causes  the  system  to  attain  position  (3) .  The 
solution  then  oscillates  between  positions  (2)  and  (3).  Whether  another 
case  using  this  configuration  will  converge  or  not  depends  on  both  the 
magnitude  of  the  current  and  the  system  weight.  The  particular  case  in 
question  probably  would  converge  in  a  more  uniform  current  field  or  if 
the  system  weight  were  to  be  increased. 

Another  non-convergent  case  involved  a  sharp  change  in  the  velocity 
profile. 


current  profile 


The  subsurface  buoy  was  initially  in  a  high  current  zone  (sketch  I);  in 
the  first  iteration  it  moved  to  a  lower  current  region  (sketch  II). 

Here  the  high  current  was  no  longer  acting  on  the  system  so  on  the  next 
iteration  the  buoy  was  moved  back  into  the  high  current  regime.  This 
case  failed  to  converge  primarily  because  the  current  shear  was  very 
near  the  initial  depth  of  the  buoy.  Thus  the  buoy  moved  in  and  out  of 
the  high  current  on  subsequent  iterations.  If  the  current  shear  had 
been  much  below  the  buoy  and  always  acted  only  on  the  cable  the  solution 
most  likely  would  have  converged. 
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The  above  cases  are  presented  simply  to  illustrate  non-convergent 
problems.  They  by  no  means  represent  all  non-convergent  cases;  however, 
they  serve  to  illustrate  what  kinds  of  situations  can  be  difficult  to 
solve.  If  a  user  encounters  a  case  that  does  not  converge,  the  discussion 
of  the  examples  here  will  help  guide  his  reasoning  in  determining  why 
non-convergence  occurs. 

In  an  attempt  to  improve  the  iteration  technique,  a  half-step 
iteration  technique  has  been  introduced.  The  process  is  implemented 
only  after  100  standard  iterations  have  occurred.  (The  100  iteration 
number  has  been  relatively  arbitrarily  judged  to  be  an  indicator  of 
oscillation  during  iteration.)  After  the  first  100  iterations  have 
occurred  the  next  100  iteration  steps  are  calculated  in  the  same  way  but 
the  allowed  displacement  is  reduced  by  one-half.  This  approach  has  not 
yet  (as  of  this  writing)  been  user  tested;  however,  initial  tests  indicate 
that  formerly  oscillating  cases  will  converge  to  a  solution.  As  part  of 
the  half-step  technique,  an  informative  printout  has  been  added  that 
states  the  number  of  iterations  required  to  achieve  convergence.  After 
200  total  iterations  have  occurred,  iteration  will  be  stopped,  a  message 
printed  and  the  next  parametric  case  will  be  analyzed.  (The  200  itera¬ 
tion  limit  is  arbitrarily  taken  as  an  indicator  that  the  solution  is  not 
converging. ) 

Convergence  Parameters 

The  variable  on  the  COMP  card  is  used  as  the  convergence  test 
parameter.  The  discussion  of  the  COMP  card  defines  a  lower  limit  for 
the  parameter  that  can  be  quite  small.  Iteration  continues  until  two 
consecutive  calculated  positions  for  each  node  on  the  structure  differ 
by  less  than  the  convergence  parameter.  Very  small  values  of  the  conver¬ 
gence  parameter  can  cause  certain  "sensitive"  cases  (such  as  some  of 
those  discussed  above)  to  appear  to  oscillate.  Adequate  solutions  and 
more  rapid  convergence  can  be  obtained  by  picking  a  convergence  parameter 
value  that  is  consistent  with  the  size  of  the  structure  being  analyzed. 

The  table  below  lists  suggested  convergence  parameter  values  that  are 
consistent  with  the  structure  size  as  determined  by  cable  length. 

Cable  Length  (ft)  Convergence  Parameter  (ft) 


<10 

<0.01 

<100 

0.01-0.1 

<1000 

0. 1-1.0 

<10,000 

1.0-10.0 

>10,000 

5.0-20.0 

These  are  only  suggested  values  that  may  be  used  as  initial  values  in 
order  to  insure  that  successful  convergence  does  occur.  For  particular 
problems  the  user  must  judge  for  himself  the  adequacy  of  the  convergence 
tolerance  chosen.  However,  the  user  must  realize  that  an  unnecessarily 
small  convergence  parameter  can  cause  a  case  to  fail  to  converge  mathe¬ 
matically  even  though  a  physically  adequate  solution  has  been  reached  in 
the  iteration  process. 


INPUT  CARDS 


Input  cards  to  DECEL1  contain  a  4-digit  integer  (-999  to  9999)  card 
number  in  columns  1-4,  a  4-character  card  type  identifier  in  columns 
5-8,  and  descriptive  properties  of  the  array  in  the  remaining  fields 
with  the  exception  of  the  following  input  cards. 

1  -  Main  Descriptive  Title  Card 

2  -  Parametric  Descriptive  Title  Card 

3  -  Continuation  Card  for  CAB 

4  -  Velocity  Field  for  Current  Option  3 

The  card  number  has  two  uses.  First,  it  is  used  as  a  convenience 
for  the  user  to  define  the  order  of  the  cards  in  the  deck,  should  the 
deck  be  dropped.  Secondly,  it  is  used  as  a  cross-check  in  changing 
parameters  in  the  Parametric  Study  Source  Deck:  when  a  parameter  on  a 
card  is  to  be  changed,  both  the  card  number  and  type  must  match  the  card 
in  the  Cable  Array  Source  Deck  or  an  error  will  occur.  In  the  Cable 
Array  Source  Deck,  duplicate  card  numbers  are  detected  as  errors;  however 
any  number  of  input  cards  may  have  the  card  number  omitted  with  no  errors 


24 


? 


CABLE  ARRAY  SOURCE  DECK  CARDS  ** 


LUN  CARD  (optional  card;  if  used  it  must  be  the  first  card  in  the  deck) 


Field 

Columns 

Format 

Contents 

Comments 

1-4 

14 

Card  number  (optional) 

-999  to  9999 

5-8 

A4 

-LUN 

Right  adjust 

| 

9-16 

18 

Unit  number  of  card  reader 

4 

17-24 

18 

Unit  number  of  line 
printer 

5 

25-32 

18 

Unit  number  of  temporary 
storage  tape 

6 

33-40 

18 

Input  option  (0  or  1) 

See  notes 

7 

41-48 

18 

Blank  if  input  option  =  0. 
Unit  number  of  source  tape 
if  input  option  =  1 . 

8 

49-56 

18 

Output  option  (0,  1  or  2) 

See  notes 

9 

57-64 

18 

Blank  if  output  option  =  0. 
Unit  number  of  output  tape 
or  card  punch  if  output 
option  =  1  or  2. 

10 

65-72 

Not  used 

See  notes 

11 

73-77 

Not  used 

12 

78-80 

Not  used 

NOTES:  1.  The  carot  symbol  (-)  is  used  to  indicate  one  blank 

column  in  the  position  shown. 

2.  The  LUN  card  is  used  to  transmit  the  logical  unit  numbers  of 
of  the  I/O  devices  and  the  I/O  options.  The  LUN  card  is 
optional.  If  it  is  omitted  then  input  option  0  and  output 
option  0  are  the  default  options. 

Two  input  options  are  available: 

0  -  The  physical  characteristics  of  the  array  are  to  be 
read  from  the  cable  array  source  deck  (see  CABLE 
ARRAY  SOURCE  DECK) 

1  -  The  physical  characteristics  of  the  array  are  to  be 
read  from  the  cable  array  source  tape  (see  CABLE 
ARRAY.  SOURCE  TAPE) 


Three  output  options  are  available: 

0  -  A  structural  output  to  the  line  printer  (see 
STRUCTURAL  OUTPUT) 

1  ”  A  device  location  output  to  tape  or  cards  (see 
DEVICE  LOCATION  OUTPUT) 

3  -  both  of  the  above 
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MAIN  DESCRIPTIVE  TITLE  CARD 


Field 

Columns 

Format 

Contents 

Comments 

1-10 

1-80 

8A10 

Main  descriptive  title 

See  notes 

NOTE:  A  single  "Main  Descriptive  Title  Card"  must  always  follow 
the  LUN  card.  If  no  title  is  desired  a  blank  card  must 
be  inserted.  This  card  does  not  adhere  to  the  convention 
requiring  a  card  number  and  card  type,  so  that  the  entire 
card  can  be  utilized  for  a  title.  If  the  LUN  card  is 
omitted  then  this  must  be  the  first  input  card. 
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NJNC  CARO 


Field 

Columns 

Format 

Contents 

Comments 

1 

1-4 

14 

Card  number  (optional) 

-999  to  9999 

2 

5-8 

A4 

NJNC 

3 

9-16 

18 

Number  of  junctions  in 
original  (unreduced) 
array 

2-44 

4 

17-24 

Not  used 

5 

25-32 

Not  used 

6 

33-40 

Not  used 

7 

41-48 

Not  used 

8 

49-56 

Not  used 

9 

57-64 

Not  used 

10 

65-72 

Not  used 

11 

73-77 

Not  used 

12 

78-80 

Not  used 

NOTE:  The  NJNC  card  is  used  to  transmit  the  number  of  junctions  in 
the  original  (unreduced  array). 
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ANC  CARO 


Field 

Columns 

Format 

Contents 

Comments 

1 

1-4 

14 

Card  number  (optional) 

-999  to  9999 

2 

5-8 

A4 

-.ANC 

Right  adjust 

3 

9-16 

18 

Junction  number  of 
anchor 

1-44 

4 

17-24 

F8.0 

X  coordinate  of  anchor 

ft 

5 

25-32 

F8.0 

Y  coordinate  of  anchor 

ft 

6 

33-40 

F8.0 

Z  coordinate  of  anchor 

ft 

7 

41-48 

Not  used 

8 

49-56 

Not  used 

9 

57-64 

Not  used 

10 

65-72 

Not  used 

11 

73-77 

Not  used 

12 

78-80 

Not  used 

NOTE:  The  ANC  cards  are  used  to  transmit  the  data  in  Table  1. 

There  must  be  one  ANC  card  for  each  anchor  in  the  array. 
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IR  CARO 


Field 

Columns 

Format 

Contents 

Comments 

1 

1-4 

14 

Card  number  (optional) 

-999  to  9999 

2 

5-8 

A4 

~IR 

Right  adjust 

3 

9-16 

18 

Junction  number 
assigned  to  a  cut  in 
the  reduced  array 

1-44 

4 

17-24 

18 

Junction  number  at 
which  cut  is  made  in 
the  original  array 

1-44 

5 

Not  used 

6 

Not  used 

7 

41-48 

Not  used 

8 

49-56 

Not  used 

9 

57-64 

Not  used 

10 

65-72 

Not  used 

11 

73-77 

Not  used 

12 

78-80 

Not  used 

NOTE:  The  IR  cards  are  used  to  transmit  the  information  contained 

in  the  table  below  the  left-hand  schematic  of  each  of  Figures 
3a-3d.  There  must  be  one  IR  card  for  each  cut  made  in  going 
from  the  original  to  the  reduced  array. 


CAB  CARD 


Field 

Columns 

Format 

Contents 

Comments 

1 

1-4 

14 

Card  number  (optional) 

-999  to  9999 

2 

5-8 

A4 

-CAB 

Right  adjust 

3 

9-16 

18 

Cable  number 

1-22 

4 

17-24 

18 

s=0  junction  number 

1-44 

5 

25-32 

18 

s=L  junction  number 
(after  required  cuts 
are  made) 

1-44 

6 

33-40 

F8.0 

Cable  weight  per 
length  in  surrounding 
fluid;  +  if  positively 
buoyant,  -  if  nega¬ 
tively  buoyant 

lb/ ft 

7 

41-48 

F8.0 

Normal  drag  coefficent 
of  cable 

8 

49-56 

F8.0 

Cable  diameter 

in. 

9 

57-64 

F8.0 

Total  cable  length 
(unstressed) 

ft 

10 

65-72 

F8.0 

Cable  rigidity,  C  (see 
notes)  if  k=l,  C=EA 

lb 

11 

73-77 

F5.0 

Exponent  in  constitu¬ 
tive  relation,  k  (see 
notes) 

>0 

12 

78-80 

13 

Number  of  straight  ele¬ 
ments  by  which  cable 
is  to  be  represented 

>0,<50 

NOTE:  The  CAB  cards  are  used  to  transmit  the  information  contained 
in  the  table  below  the  right-hand  schematic  in  Figures  3a-3d 
(fields  3  to  5),  the  physical  characteristics  of  the  cables 
in  the  array  (fields  6  to  11),  and  the  fineness  by  which  the 
cables  in  the  array  are  to  be  modeled  (field  12).  There 
must  be  one  CAB  card  and  one  continuation  CAB  card  for  each 
cable  in  the  array.  The  constitutive  relation  for  a  cable 
can,  in  general,  be  written  (6)  as  e=(T/C;  where  e  is  the 
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NOTES  FOR  CAB  CARD  (continued) 


strain  and  T  the  tension.  C  and  k  are  constants  which 
depend  on  the  cable  material  and  construction.  Fields  10 
and  11  transmit  the  values  of  C  and  k,  respectively.  An 
inextensible  cable  is  transmitted  to  DECEL1  by  leaving 
fields  10  and  11  blank. 
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CONTINUATION  CAB  CARD 


Field 

Columns 

Format 

Contents 

Comments 

1 

Tangential  drag  coef¬ 
ficient  of  cable 

_ | 

See  notes 

2 

9-80 

9F8.0 

Not  used 

NOTE:  Each  CAB  card  must  be  followed  by  an  additional  card.  If 

no  tangential  drag  is  desired  a  blank  card  must  be  inserted. 
This  card  does  not  adhere  to  the  convention  requiring  a  card 
number  and  card  type  since  it  is  in  fact  a  continuation  of 
the  tangential  drag  coefficient.  Default  value  is  zero. 
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DCAB  CARD 


Field 

Columns 

Format 

Contents 

Comments 

1 

1-4 

14 

Card  number  (optional) 

-999  to  9999 

2 

5-8 

A4 

DCAB 

i 

3 

9-16 

18 

Number  of  cable  to 
which  device  is 
attached 

1-22  i 

i 

4 

17-24 

18 

Device  type  (1,2) 

] 

See  notes 

-  -i 

5 

25-32 

18 

Print,  flag  for  cable 
devices.  If  flag  =  0, 
print  device  charac¬ 
teristics  and  location. 

If  flag  f  0,  don't 
print. 

; 

i 

6 

33-40 

F8.0 

Device  weight  in  sur¬ 
rounding  fluid;  +  if 
positively  buoyant,  - 
if  negatively  buoyant 

ib 

(See  note  3) 

7 

41-48 

F8.0 

Normal  drag  coefficient 

8 

49-56 

F8.0 

Diameter  on  which  drag 
coefficient  is  based 
(Type  1) 

in. 

Frontal  area  on  which 
drag  coefficient  is 
based  (Type  2). 

ft2 

9 

F8.0 

Device  length  (Type  1) 

ft 

10 

65.-72 

F8.0 

Unstressed  distance  of 
device  from  s=0  junc¬ 
tion  of  the  cable 

ft 

11 

73-77 

F5.0 

Tangential  drag  coef¬ 
ficient 

12 

78-80 

Not  used 

NOTE:  1.  The  DCAB  cards  are  used  to  transmit  the  physical  charac¬ 

teristics  of  the  discrete  devices  (buoyancy  elements, 
current  meters,  etc.)  attached  to  the  cables  in  the 
airay.  There  must  be  one  DCAB  card  for  each  such  device 


NOTES  FOR  DCAB  CARD  (continued) 


in  the  array.  Two  types  of  devices  are  permitted:  (a) 
in-line  (Type  1)  any  elongated  device  (cylinder,  ellipse, 
etc.)  attached  so  that  its  longitudinal  axis  is  aligned 
with  the  cable  axis;  and  (b)  free  devices  (Type  2). 

2.  Dummy  DCAB  devices  with  no  weight  and  zero  drag  can  be 
used  as  conveniences  to  print  the  location  of  the  device 
in  the  cable  array's  deformed  state.  With  no  DCAB 
devices,  no  information  about  the  cable  shape  between 
junctions  is  printed. 

3.  Include  cable  weight  in  water  if  device  covers  segment 
of  cable. 

Each  device  is  automatically  assigned  a  unique  number  based  on 
its  location  on  the  structure.  Low  numbered  devices  are  located 
on  low  numbered  cables.  On  each  cable  the  lowest  number  device 
is  nearest  the  s=o  end  of  the  cable. 
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DJNC  CARD 


Field 

Columns 

Format 

Contents 

Comments 

1 

1-4 

14 

Card  number  (optional) 

-999  to  9999 

2 

5-8 

A4 

DJNC 

3 

9-16 

18 

Number  of  junction  to 
which  device  is 
attached 

1-44 

4 

17-24 

18 

Not  used 

5 

25-32 

18 

Not  used 

6 

33-40 

F8.0 

Device  weight  in  sur¬ 
rounding  fluid;  +  if 
positively  buoyant,  - 
if  negatively  buoyant 

lb 

7 

41-48 

F8.0 

Device  drag  coefficient 

8 

49-56 

F8.0 

Frontal  area  on  which 
drag  coefficient  is  based 

ft2 

9 

57-64 

Not  used 

10 

65-72 

Not  used 

11 

73-77 

Not  used 

12 

78-80 

Not  used 

NOTE:  The  DJNC  cards  are  used  to  transmit  the  physical  character¬ 
istics  of  the  discrete  devices  attached  to  the  junctions  in 
the  array.  There  must  be  one  DJNC  card  for  each  such  device 
in  the  array. 

Since  an  in-line  device  cannot  physically  exist  at  a  junc¬ 
tion  (cable  termination),  only  free  devices  are  permitted 
at  a  junction. 

DJNC  devices  are  indexed  automatically  for  counting  pur¬ 
poses  only;  they  are  counted  in  the  total  number  of  indexed 
devices . 
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DEN  CARD 


Field 

Columns 

Format 

1  Contents 

Comments 

9-16 

17-24 

8 

49-56 

9 

57-64 

10 

65-72 

11 

73-77 

12 

78-80 

Card  number  (optional) 


A4  .DEN 


which  array  is 
suspended 


Not  used 


Not  used 


Not  used 


Not  used 


Not  used 


Not  used 


Not  used 


Not  used 


Not  used 


-999  to  9999 


Right  adjust 


lb  s2  ft"4 
(1.99  for 
seawater) 


NOTE:  The  DEN  card  is  used  to  transmit  the  density  of  the  fluid  in 
which  the  array  is  suspended. 


EOD  CARD 


Field 

Columns 

Format 

Contents 

Comments 

1-4  14  Card  number  (optional) 


5-8  A4  -EOD 


9-16  Not  used 


17-24  Not  used 


25-32  Not  used 


33-40  Not  used 


41-48  Not  used 


49-56  Not  used 


57-64  Not  used 


65-72  Not  used 


73-77  Not  used 


78-80  Not  used 


-999  to  9999 


Right  adjust 


The  EOD  card  is  used  to  specify  the  end  of  data  transmission. 


NDAT  CARD 


Field 

Columns 

Format 

Contents 

Comments 

1 

— 

1-4 

14 

Card  number  (optional) 

-999  to  9999 

2 

5-8 

A4 

NDAT 

3 

9-16 

18 

Current  field  option 

See  note  1 

4 

17-24 

18 

Number  of  stations  at 
which  current  data  is 
taken  for  current  field 
option  3.  Otherwise 
not  used. 

See  note  2 

5 

25-32 

18 

Option  for  defining  a 
particular  parametric 
study  as  the  new  no 
current  reference  con¬ 
figuration.  Option  = 

0,  no  new  ref.  config. 
Option  t  0,  new  ref. 
config. 

6 

33-40 

F8.0 

PHI 

degrees 

See  note  3 

7 

41-48 

F8.0 

Device  location  punch 
output  option;  option  / 

0,  punch  device  loca¬ 
tions  on  cards 

See  note  4 

8 

49-56 

F8.0 

Weight  of  device  whose 
location  is  to  be 
punched 

lb 

9 

57-64 

Not  used 

10 

65-72 

F8.0 

Velocity  units  input 
option 

See  note  5 

11 

73-77 

F5.0 

Number  of  steps  used  to 
apply  full  current 

See  notes  6 
and  7  (default 
value  =  1.0) 

12 

78-80 

Not  used 

SEE  NEXT  PAGE  FOR  NOTES 
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NOTE: 

1.  The  NDAT  card  is  used  to  specify  that  new,  modified,  or  addi¬ 
tional  data  follow  and  to  transmit  the  current  field  option. 

Three  current  field  options  are  available. 

0  -  No  current. 

1  -  Current  field  is  unidirectional  and  horizontal. 

(There  is  no  vertical  component.)  The  current  magni¬ 
tude  can  vary  with  depth.  At  any  point  in  the  plane 
Z-constant,  the  velocity  is  invariant.  The  no-current 
configuration  is  also  calculated  as  part  of  this 
option. 

2  -  Same  as  current  field  option  1  except  that  the  current 

can  have  a  directionality  that  depends  on  depth.  That 
is,  this  option  relaxes  the  unidirectional  constraint 
on  the  velocity  field.  At  any  point  in  the  plane 
Z-constant,  the  velocity  is  invariant.  An  example  of 
a  possible  current  candidate  requiring  this  option 
would  be  a  helical  current  where  the  axis  of  the  helix 
is  aligned  with  gravity.  The  no-current  configuration 
is  also  calculated  as  part  of  this  option. 

3  -  Interpolation/extrapolation  scheme  for  field  current 

(magnitude  and  direction)  data.  The  current  is 
assumed  to  have  no  vertical  component.  Two,  three, 
or  four  stations  can  be  selected  in  the  horizontal 
plane  and  at  each  of  these  stations,  current  (direc¬ 
tion  and  magnitude)  data  can  be  obtained  at  up  to  25 
elevations.  This  option  affords  current  variation 
in  horizontal  planes.  The  no-current  configuration 
is  also  calculated  as  part  of  this  option. 

2.  For  current  option  3,  the  value  in  field  4  on  the  NDAT  card  must 
appear  on  all  successive  NDAT  cards  even  if  the  velocity  field 
is  not  to  be  varied  in  any  of  the  parametric  studies. 

3.  The  angle  0  specifies  the  rotation  of  the  X-axis  from  the  N-axis. 

<|>  is  measured  positive  in  the  clockwise  sense  and  the  N-axis  is 
0=0.  That  is,  when  0=0  then  the  (N,W)  coordinate  system  is  coin¬ 
cident  with  the  (X,Y)  coordinate  system  (see  Figure  4).  The  0 
value  should  appear  on  each  NDAT  card. 

4.  It  is  sometimes  desirable  to  have  the  locations  of  particular 
devices  punched  on  cards  for  input  to  other  programs.  For  exam¬ 
ple,  both  hydrophones  and  buoys  may  be  distributed  as  DCAB  devices 
along  the  structure,  but  only  the  location  of  the  hydrophones  is 
important.  By  specifying  field  7  f  0  and  field  8  =  hydrophone 
weight,  DECEL1  will  punch  for  each  hydrophone:  (a)  device  index, 
(b)  cable  number,  (c)  distance  from  the  s  =  0  end  of  the  cable, 
and  (d)  x,y,z  coordinates  of  the  device.  Preceding  this  punched 
output  is  a  card  totalling  the  number  of  cards  in  the  punched 
output. 
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NOTES  FOR  NDAT  CARD  (continued) 

5.  Velocity  units  input  option: 

0  -  input  magnitude  of  velocity  in  knots  (default) 

1  -  input  magnitude  of  velocity  in  cm/sec 

2  -  input  magnitude  of  velocity  in  ft/sec 

6.  This  is  a  coded  value:  the  integer  portion  represents  the  num¬ 

ber  of  steps  used  to  apply  the  full  current;  the  fractional 
part  controls  printed  output:  a  non-zero  fractional  part  (i.e., 
10.1)  causes  the  shape  calculated  for  each  increment  of  the 
current  to  be  printed.  For  example:  the  value  10.1  will  cause 
the  current  to  be  applied  in  ten  equal  increments  and  the  struc¬ 
ture  shape  will  be  printed  for  each  increment.  A  10.0  value 
will  apply  the  current  in  ten  increments  but  print  only  the 
shape  with  the  full  current  value  applied.  If  plotting  is 
requested,  only  the  shape  with  the  full  current  applied  is 
plotted. 

7.  This  input  option  is  intended  to  aid  in  obtaining  solutions 
where  high  currents  cause  SLOW/NO  CONVERGENCE  messages  to  be 
printed. 
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PARAMETRIC  DESCRIPTIVE  TITLE  CARD 


NOTE:  Descriptive  Title  Card  must  always  follow  an  NDAT  card.  If 
no  title  is  desired,  then  a  blank  card  must  be  inserted. 

This  card  does  not  adhere  to  the  convention  requiring  a  card 
number  and  card  type,  so  that  the  entire  card  can  be  utilized 
for  a  title. 
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COMP  CARD 


Field 

Columns 

Format 

Contents 

Comments 

1 

1-4 

mm 

Card  number  (optional) 

-999  to  9999 

2 

5-8 

A4 

COMP 

3 

9-16 

F8.0 

Accuracy  required  in 
array  equilibrium  cal¬ 
culations  (COMP  value) 

ft 

(0.01  typ) 

4 

17-24 

F8.0 

Nodal  position  itera¬ 
tion  limit 

5 

25-32 

F8.0 

Imaginary  reaction 
iteration  limit 

See  note  2 
(1000  default 
value) 

6 

33-40 

Not  used 

1 

7 

41-48 

Not  used 

8 

49-56 

Not  used 

9 

57-64 

Not  used 

10 

65-72 

Not  used 

11 

73-77 

Not  used 

12 

78-80 

Not  used 

NOTES: 


1.  The  COMP  card  is  used  to  transmit  the  accuracy  requirement 
imposed  on  the  array  equilibrium  calculations.  This  accu¬ 
racy  requirement,  specified  in  field  3,  insures  that  the 
calculated  coordinates  of  every  point  in  the  array  are 
within  ±  field  3  of  their  exact  values. 

The  accuracy  obtainable  is  limited  by  the  significant  figure 
capacity  of  the  computer  being  used  and  by  the  largest 
linear  dimension  in  the  array.  Let  the  number  of  signifi¬ 
cant  figures  carried  in  single  precision  be  n,  and  let  the 
characteristic  of  the  common  logarithm  of  the  largest  linear 
dimension  be  m.  Then,  the+galue  of  field  3  is  usually 
bounded  by  field  3  >  10  .  (For  example,  suppose  n  =  8 

and  the  largest  dimension  is  25,000  ft.  Then,  field  3  >  0.1 
ft.)  Occasionally,  a  larger  minimum  value  must  be  used.  A 
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NOTES  FOR  COMP  CARD  (continued) 

COMP  card  must  appear  after  the  first  NDAT  card.  See  PARA 
METRIC  STUDY  SOURCE  DECKS.  Accuracy  requirements  can  be 
changed  in  subsequent  Parametric  Study  Decks  by  using  new 
COMP  cards. 

2.  See  the  descriptions  in  the  text  on  pages  18  and  19. 
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VEL  CARD 


Field  Columns  I  Format 


Contents 


Comments 


Z  coordinate  at  which 
current  velocity  is 
specified 


Magnitude  of  current 
at  Z  coordinate  speci 
fied  in  field  3 


Direction  of  current 
at  Z  coordinate  (cur¬ 
rent  option  2  only) 


Not  used 


-999  to  9999 

Right  adjust 

ft 

knots ; 
cm/sec; 
ft/sec  (con¬ 
sistent  with 

NDAT  card) 

degrees 

Not  used 


Not  used 


Not  used 


Not  used 


NOTE:  The  VEL  cards  are  used  to  transmit  the  data  in  Table  2. 

There  must  be  one  VEL  card  for  each  value  of  Z  at  which  the 
current  velocity  is  specified. 

Up  to  25  VEL  cards  are  permitted.  At  least  one  of  the  Z 
coordinates  transmitted  by  the  VEL  cards  must  be  less  than 
or  equal  to  the  minimum  Z  coordinate  transmitted  by  the  ANC 
cards.  For  further  information  see  CURRENT  FIELD  INPUT 
OPTION  1. 

The  VEL  cards  must  not  appear  after  an  NDAT  card  specifying 
a  current  field  option  =  0.  All  VEL  cards  required  to 
transmit  the  current  profile  must  appear  after  the  first 
NDAT  card  specifying  a  current  field  option  =  1.  See  PARA¬ 
METRIC  STUDY  SOURCE  DECKS. 

A  given  current  profile  will  be  applied  in  all  subsequent 
parametric  cases  until  a  new  field  is  defined. 
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NOTES  FROM  VEL  CARD  (continued) 


Current  Option  1  -  Only  the  magnitude  of  the  current  is 

input  (field  4).  The  direction  is  varied  by  the  ANG 
card . 

Current  Option  2  -  Both  magnitude  and  direction  of  current 
are  input  fields  4  and  5,  respectively.  The  direction 
of  the  field  also  may  be  modified  by  the  ANG  card. 
Current  direction  is  to  be  specified  from  the  magnetic 
north  axis.  Clockwise  rotation  is  positive. 

Current  Option  3  -  Fields  3,  4  and  5  are  ignored.  The  velo¬ 
city  information  is  read  from  the  velocity  continuation 
group  immediately  following  the  VEL  card. 

The  following  is  a  brief  description  of  the  velocity 
cards  under  current  option  3. 

After  a  VEL  card  is  encountered  there  will  be  2,  3  or  4 
sets  of  velocity  data  depending  on  the  number  of  sta¬ 
tions  at  which  velocity  profiles  were  measured. 

Each  set  of  data  for  a  station  will  contain  one  Station 
Location  Card  and  N  Velocity  Definition  Cards.  N  is 
the  number  of  Z-coordinates  at  which  velocity  is  mea¬ 
sured  for  the  particular  station.  The  maximum  value  of 
N  is  25. 

The  2,  3  or  4  sets  of  data  are  stacked  as  indicated  in 
Figure  7. 

The  Station  Location  Cards  and  Velocity  Definition 
Cards  are  defined  on  the  next  two  pages. 
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To  further  clarify  the  velocity  option  3  data  requirements  we  present  another 
example  using  3  stations  of  current  data.  The  data  deck  will  be  arranged  as 
follows: 


I 


Figure  7.  Arrangement  of  data  sets  for  velocity  option  3. 


NOTES  FROM  VEL  CARD  (OPTION  3  ONLY)  (continued) 


Field 

Columns 

Format 

Contents 

Comments 

1 

- 

1-5 

15 

Station  number  at  which 
velocity  profile  mea¬ 
sured 

2 

6-10 

15 

Number  of  Z-coordi- 
nates  at  which  velo¬ 
city  is  measured  for 
this  station  (field  1) 

3 

11-20 

F10.0 

North  coordinate  of 
this  station* 

ft 

(see  notes) 

4 

21-30 

F10.0 

West  coordinate  of 
this  station** 

ft 

(see  notes) 

5 

31-80 

Not  used 

♦Coordinate  is  positive  for  north,  negative  for  south. 
♦♦Coordinate  is  positive  for  west,  negative  for  east. 


50 


NOTES  FROM  VEL  CARD  (OPTION  3  ONLY)  (continued) 


Field 

Columns 

Format 

Contents 

Comments 

1 

1-5 

15 

Station  number 

2 

6-10 

Not  used 

3 

11-20 

F10.0 

Z-coordinate  of  velo¬ 
city 

ft 

4 

21-30 

F10.0 

Magnitude  of  velocity 
at  Z-coordinate  of  this 
station  (field  1)* 

See  notes 

5 

31-40 

F10.0 

Direction  of  current 
at  Z-coordinate  of 
this  station  (field 

1)** 

degrees 
(see  notes) 

6 

41-80 

Not  used 

NOTE:  There  must  be  as  many  of  this  card  as  there  are  current  read¬ 
ings  at  the  station  (=  field  2  of  preceding  card) . 


*The  magnitude  of  velocity  may  be  input  in  three  different  units 
according  to  the  option  definition  on  the  NDAT  card. 

Units  Option  0  -  knots  (default) 

Units  Option  1  -  cm/sec 

Units  Option  2  -  ft/sec  (see  NDAT  card) 


**The  direction  of  current  is  positive  in  the  clockwise  sense  from 
the  magnetic  north  axis  (degrees). 
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ANG  CARD 


Field 

Columns 

Format 

Contents 

Comments 

1 

1-4 

14 

Card  number  (optional) 

-999  to  9999 

2 

5-8 

A4 

-ANG 

Right  adjust 

3 

9-16 

F8.0 

Initial  current  angle 

degrees 

4 

17-24 

F8.0 

Increment  in  current 
angle 

deg,  >  0 

5 

25-32 

F8.0 

Final  current  angle 

deg, 

>  field  3 

6 

33-40 

Not  used 

7 

41-43 

Not  used 

8 

49-56 

Not  used 

9 

57-64 

Not  used 

10 

65-72 

Not  used 

11 

73-77 

Not  used 

12 

78-80 

Not  used 

NOTE:  The  ANG  card  is  used  to  transmit  instructions  for  changing 
the  angular  direction  of  the  current  with  respect  to  the 
N-axis  for  current  options  1  and  2  only.  DECEL1  calculates 
the  array  response  to  the  specified  current  profile  from  the 
initial  to  the  final  current  angles  in  increments  trans¬ 
mitted  by  field  4.  An  ANG  card  must  not  appear  after  an 
NDAT  card  specifying  a  current  field  option  =  0.  One  ANG 
card  must  appear  after  the  first  NDAT  card  specifying  a 
current  field  option  =  1,  2,  or  3.  See  PARAMETRIC  STUDY 
SOURCE  DECKS. 

If  no  parametric  range  of  current  direction  variation  is 
required,  then  Fields  3  and  5  should  have  identical  values 
and  field  4  should  have  a  non-zero  positive  value. 

When  exercising  the  angular  rotation  option,  it  is  important 
to  recognize  the  difference  between  current  options  1  and  2. 
For  current  option  1,  the  input  value  in  field  3  of  the  ANG 
card  establishes  the  directionality  of  the  entire  flow  field 
This  is  not  the  case  for  current  option  2  since  current 
directionality  is  established  and  specified  as  input  on  the 
VEL  card. 
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NOTES  FROM  ANG  CARD  (continued) 


For  current  option  2  the  field  3  input  value  on  the  ANG 
card  should  be  set  equal  to  $  (the  angle  between  the  N- 
and  X-axes)  plus  the  initial  angle  of  interest.  The  value 
in  field  5  of  the  ANG  card  should  be  set  equal  to  $  +  0 
where  p  is  the  total  angle  through  which  the  current  is  to 
be  rotated.  The  value  of  field  4  is  Ap,  the  increment. 
When  current  option  2  is  selected  and  no  rotation  of  the 
current  field  is  desired,  the  fields  3  and  5  should  be  set 
equal  to  <}>  and  any  positive  value  set  in  field  4. 


PPLT  CARD 


Field 

Columns 

Format 

Contents 

Comments 

1 

1-4 

14 

Card  number  (optional) 

-999  to  9999 

2 

5-8 

A4 

PPLT 

3 

9-16 

F8.0 

Plotting  option: 

0  -  Undeformed  only 

1  -  Deformed  only 

2  -  Both 

Undeformed 
plotted  as  a 
dotted  line. 
Deformed 
plotted  as  a 
solid  line. 

4 

17-24 

F8.0 

Height  of  plot,  y 

in. 

(default  10) 

5 

25-32 

F8.0 

View  angle,  x 
(see  note) 

degrees 
(default  30) 

6 

33-40 

F8.0 

View  angle,  y 
(see  note) 

degrees 
(default  120) 

7 

41-48 

F8.0 

View  angle,  Z 
(see  note) 

_ 

degrees 
(default  90) 

NOTE:  For  a  plan  or  elevation  view  (not  a  perspective  view)  one 
view  angle  must  be  361°.  This  is  a  code  indicating  which 
axis  is  perpendicular  to  the  plot.  For  example,  a  plan  view 
is  specified  by: 

View  Angle,  x  =  0.001  (0  gives  the  30°  default  value) 
View  Angle,  y  =  90° 

View  Angle,  z  =  361° 

The  x,  y  and  z  view  angles  are  shown  below  for  the  default 
configuration. 


z 


9^  -  View  angle  for  x  axis.  30  degrees. 

8^  -  View  angle  for  y  axis.  120  degrees. 

0z  -  View  angle  for  z  axis.  90  degrees. 

Examples  perspective  plots  are  shown  in  Figures  8-10. 
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CPLT  CARD 


Field 

Columns 

Format 

Contents 

Comments 

1 

1-4 

14 

Card  number  (optional) 

-999  to  9999 

2 

5-8 

A4 

CPLT 

3 

9-16 

F8.0 

ZUP-Max.  depth  plane 
level 

ft 

1  4 

17-24 

J 

F8.0 

ZDN-Min.  depth  plane 
level 

ft 

5 

25-32 

F8.0 

DZ  -Delta  depth 

ft 

6 

33-40 

F8.0 

YIN-Height  in  Y  direc¬ 
tion  of  a  single  plane 

in. 

7 

41-48 

F8.0 

XMIN  | 

ft 

(If  blank 
default 
values 
selected) 

8 

49-56 

F8.0 

/  Minimum  &  maxi- 
XMAX  f  ,  . 

1  mum  coord,  of 

9 

57-64 

F8.0 

/  area  to  be  cov- 
YnlN  l  ,  ,  ,  . 

\  ered  by  plot 

10 

65-72 

F8.0 

. 1 

YMAX  ) 

11 

73-77 

F5.0 

ANG-View  angle  Y 

degrees 
default  =  90 

12 

78-80 

13 

NY  -Number  of  mesh 

points  in  Y  direc¬ 
tion,  including 
boundary 

default  =  6 

NOTE:  Current  plots  may  be  depicted  at  one  or  more  depths  by  vary¬ 
ing  input  parameters  ZUP,  ZDN,  DZ.  The  vertical  height  for 
each  plane  of  the  plot  is  input  as  YIN,  and  the  correspond¬ 
ing  width  is  calculated  by  the  program.  The  product  of  the 
number  of  planes  and  the  vertical  height  of  a  single  plane 
cannot  exceed  10  inches . 

The  plots  may  vary  in  perspective  view  angle  y  over  the 
range  0<ANG<180°,  where  ANG  =  90°  is  the  plan  view.  For 
purposes  of  plotting  the  current  field,  each  plane  is 
assumed  to  have  z  coordinate  zero  and  view  angle  x  equal 
to  zero. 


NOTES  FROM  CPLT  CARD  (continued) 


The  mesh  for  plotting  current  is  controlled  by  input  NY. 

NY  is  the  number  of  mesh  points  in  the  y  direction,  for  a 
plane,  including  the  boundary.  The  corresponding  number  of 
mesh  points,  for  the  x  direction,  is  calculated  by  the  pro¬ 
gram.  (NY  =  2  +  number  of  current  arrows  encountered  when 
moving  from  YMIN  to  YMAX.) 

The  area  covered  by  the  current  field  plot  may  be  selected 
in  two  ways.  If  the  values  XMIN,  XMAX,  YMIN,  YMAX  are  left 
blank,  the  program  determines  the  area  to  be  plotted.  The 
area  is  based  on  the  maximum  and  minimum  anchor  point  coor¬ 
dinates.  The  second  method  is  to  define  XMIN,  XMAX,  YMIN, 
YMAX  in  terms  of  the  x,  y  coordinate  system.  If  there  is 
only  one  anchor  point  XMIN,  XMAX,  YMIN,  YMAX  must  be  defined 
by  the  user.  A  star  will  be  plotted  at  each  anchor  point 
within  the  defined  area. 

Examples  of  current  field  plots  are  shown  in  Figures  11  and 

12. 
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EOD  CARD 


Field 

Columns 

Format 

Contents 

Comments 

1 

1-4 

14 

Card  number  (optional) 

-999  to  9999 

2 

5-8 

A4 

^EOD 

Right  adjust 

3 

9-16 

Not  used 

4 

17-24 

Not  used 

5 

25-32 

Not  used 

6 

33-40 

Not  used 

7 

41-48 

Not  used 

8 

49-56 

Not  used 

9 

57-64 

Not  used 

10 

65-72 

Not  used 

11 

73-77 

Not  used 

12 

78-80 

Not  used 

_ _ 

NOTE:  The  EOD  card  is  used  to  specify  the  end  of  data  transmission 
or  end  of  a  parametric  case. 
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Cable  Array  Source  Deck 


The  cable  array  source  deck  contains  all  of  the  cards  required  to 
describe  the  physical  characteristics  of  the  cable  array  under  consider¬ 
ation.  The  last  card  in  the  cable  array  source  deck  must  be  an  EOD  card 
to  signify  the  end  of  data  transmission.  The  form  of  the  cable  array 
source  deck  is  given  below. 


Cable  Array  Source  Deck 


Card  Type* 


Comments 


NJNC 

ANC 

IR 

CAB 

DCAB 

DJNC 

DEN 

EOD 


One  card 

One  for  each  array  anchor 

One  for  each  cut  going  from  original  to  reduced 
array 

One  (plus  one  continuation  card)  for  each  array  cable 
One  for  each  discrete  device  attached  to  a  cable 
One  for  each  discrete  device  attached  at  a  junction 
One  card 

Must  be  last  card  in  source  deck 


Cable  Array  Source  Tape 

The  cable  array  source  tape  is  an  alternate  means  of  transmitting 
to  DECEL1  the  physical  characteristics  of  the  cable  array  under  considera¬ 
tion.  A  recommended  program  for  generating  the  cable  array  source  tape 
from  the  cable  array  source  deck  is  given  below. 


Cable  Array  Source  Tape  Generation 


PR0GRAM  STAPE 
DIMENSI0N  DATA(IO) 

DATA(IREAD  =  Unit  number  of  card  reader) 
DATA(ITAPE  =  Unit  number  of  source  tape) 

1  READ(IREAD, 10) (DATA(I) ,1=1 , 10) ,EX,NSEG 
WRITE ( ITAPE ,11) (DATA ( I ) , 1= 1 , 1 0 ) , EX , NSEG 
IF (DATA ( 2 ) . NE . 4HDCAB ) G0  T0  3 
READ ( IRE AD , 20 )TANDRG 
WRITE ( ITAPE ,21) TANDRG 

(continued) 


*These  cards  may  be  arbitrarily  ordered  except  for  the  EOD  card.  It  is 
strongly  recommended  that  the  cards  in  the  cable  array  source  deck  be 
given  unique  card  numbers  (Field  1  of  each  card).  Note  that  LUN,  NDAT, 
COMP,  VEL,  ANG,  and  EOP  cards  are  not  permitted  in  the  cable  array 
source  deck. 
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Cable  Array  Source  Tape  Generation  (Continued) 


3  CONTINUE 

IF(DATA(2) .EQ.4H  EOD)G0  T0  2 
G0  T0  1 

10  F0RMAT (F4 . 0 , A4 , 8F8 . 0 , F5 . 0 , I 3) 

11  F0RMAT(F4.O,A4,8E15.8,/,E12.5,I3) 

20  F0RMAT(F8.O) 

21  F0RMAT(E15.8) 

2  REWIND  TAPE 

END 


Parametric  Study  Source  Decks 

The  parametric  study  source  decks  are  used  to  transmit  to  DECEL1 
accuracy  requirements,  current  fields,  and  changes  in  the  physical 
properties  of  the  cable  array  under  consideration.  Each  parametric 
study  source  deck  must  begin  with  an  NDAT  card  and  a  Parametric  Descrip¬ 
tive  Title  card  and  end  with  an  EOD  card.  The  form  of  the  parametric 
study  source  decks  is  given  below. 
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Parametric  Study  Source  Decks 


Card  Type* 


Comments 


NDAT 


Parametric 
Title  Card 

COMP 


VEL 


ANG 


ANC 

CAB 

DCAB 

DJNC 

EOD 


Must  be  first  card  in  parametric  study  source  deck. 
Field  3  (current  option)  can  change  as  0  1,  0  2, 
or  0  3.  The  changes  12,  13,  and  2  3,  even  with 
intermediate  zeros,  are  not  permitted. 

NDAT  card  must  be  followed  immediately  by  para¬ 
metric  title  card. 

Must  appear  after  first  NDAT  card.  The  accuracy 
requirement  transmitted  is  retained  until  the 
appearance  of  a  COMP  card  in  another  parametric 
study  source  deck. 

Must  not  appear  after  NDAT  card  containing  current 
option  0.  All  VEL  cards  required  to  transmit  a 
current  profile  must  appear  after  the  first  NDAT 
card  containing  a  current  option  1  or  2  or  3.  The 
current  profile  transmitted  is  retained  until  the 
appearance  of  a  VEL  card  in  another  parametric 
study  source  deck.  The  appearance  of  the  first 
VEL  card  in  a  parametric  study  source  deck  zeros 
the  entire  current  profile.  Thus,  to  change  from 
one  current  profile  to  another,  all  VEL  cards 
required  to  specify  the  new  profile  must  appear  in 
the  appropriate  parametric  study  source  deck. 

Must  not  appear  after  NDAT  card  containing  current 
option  0.  One  ANG  card  must  appear  after  the 
first  NDAT  card  containing  current  option  1,  2  or 
3.  The  current  angles  transmitted  are  retained 
until  the  appearance  of  an  ANG  card  in  another 
parametric  study  source  deck. 

Used  to  change  anchor  data.  See  Note. 

Used  to  change  cable  data.  See  Note. 

Used  to  change  discrete  device  on  cable  data.  See 
Note. 

Used  to  change  discrete  device  at  junction  data. 

See  Note. 

Must  be  last  card  in  parametric  study  source  deck. 


*These  cards  may  be  arbitrarily  ordered  except  for  the  NDAT,  Parametric 
Title,  and  EOD  cards.  Note  that  LUN,  JNC,  IR,  DEN,  and  EOP  cards  are 
not  permitted  in  the  parametric  study  source  decks. 


66 


NOTE:  The  array  design  changes  which  are  permitted  are  those 

changing  the  physical  data  of  the  array  but  not  the  over¬ 
all  geometric  layout  of,  or  the  number  of  discrete 
devices  on,  the  array.  These  changes  are  keyed  by  match¬ 
ing  the  card  number  and  type  appearing  in  a  parametric 
study  source  deck  to  the  card  number  and  type  appearing 
in  changes  in  the  array  physical  data  is  given  below. 


Summary  of  Changes  in  the  Array  Physical  Data  Permitted  (P) 
and  Not  Permitted  (NP) 


Card  Type 


1 

2 

3 

ANC 

NP 

NP 

NP 

CAB 

NP 

NP 

NP 

CONTINUATION  CAB 

P 

DCAB 

NP 

NP 

NP 

DJNC 

NP 

NP 

NP 

Field  Number 

4  5  6  7  8  9  10  11  12 

P  P  P  Not  used 

NPNPPPPPP  P  P 

1*2  NP  P  P  P  P  P  P  P 

NP  NP  P  P  P  Not  used 


OVERALL  INPUT  DECK 

The  overall  input  deck  consists  of  a  LUN  card  (optional)  specifying 
the  I/O  options  and  the  logical  unit  numbers  of  the  required  I/O  devices, 
followed  by  the  cable  array  source  deck  (or  tape) ,  followed  by  any 
number  of  parametric  study  source  decks,  and  ended  by  an  EOP  card  signi¬ 
fying  the  end  of  the  problem.  The  form  of  the  overall  input  deck  is 
given  below. 

Overall  Input  Deck 

•  LUN  Card  (optional) 

Cable  Array  Source  Deck  (or  Tape) 

Parametric  Study  Source  Decks 

•  EOP  Card 
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EXAMPLE 


1.  Compilation  of  cable  array  source  deck  errors  only. 

1000  ..LUN.. ...... ....60  ....... ....61.... -24  (or  omit  this  card) 

*****Main  Title  Card***** 

Cable  Array  Source  Deck  (or  Tape) 

100UE0P 

2.  Deflections  of  the  source  deck  cable  array  due  to  a  cur¬ 
rent  profile. 

(LUN  card  omitted) 

*****Main  Title  Card***** 

Cable  Array  Source  Deck  (or  Tape) 

1001NDAT . 1 

*****Parametric  Study  1  Title  Card***** 

1002C0MP 

1003-.VEL 

1010..  VEL 
1011. ANG 

1012.. EOD 
1013-.E0P 

3.  Deflections  of  the  source  deck  cable  array  due  to  a  cur¬ 
rent  profile  and  effects  of  buoyancy  changes  on  these 
deflections . 

(LUN  card  omitted) 

*****MAIN  TITLE  CARD***** 

Cable  Array  Source  Deck  (or  Tape)  containing  a  card: 

A  ..9DJNC.  10000  .  A  A  A  A  1 . 95aAAA  35.3 

1001NDAT . 1 

*****Parametric  (see  note)  Study  1  Title  Card***** 

1002C0MP 

1003-.VEL 

1010..  VEL 

1011..  ANG 

1012-.E0D 
1013NDAT . 1 

*****Parametric  Study  2  Title  Card***** 

A  A  a9DJNCa aaaaaa3aaaaaaaaaaaAAaaaaa12500.aa.a1 . 95aaAA  35  •  3 

1014.E0D 
1015 aEOP 

NOTE:  In  each  case  unique  card  numbers  were  assigned,  other¬ 
wise  an  error  would  have  resulted.  However,  any  or  all 
card  numbers  could  have  been  omitted  and  no  error  would 
result. 
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ERROR  MESSAGES 


DECELl  contains  a  series  of  internal  error  checks  to  insure  that  the 
original  cable  array  is  properly  reduced  to  a  statically  determinate 
array;  the  input  data  are  consistent;  and  the  deck  is  properly  structured. 

If  errors  are  found,  the  entire  list  of  input  cards  is  printed;  and 
cards  with  errors  are  identified  by  an  error  code  number.  The  text  of 
the  coded  error  message  is  printed  after  the  card  listing.  All  cards 
are  scanned  for  errors;  however,  only  the  first  error  on  a  card  is 
detected.  A  card  with  a  DECELl  detected  error  should  be  scanned  visually 
to  check  for  other  errors. 


DEFINITION  OF  ERRORS 
Type  0  - 

A  type  0  error  indicates  that  the  card  type  identifier 
(Field  2)  is  not  recognizable. 

Type  1A  -  LUN 

Field  6  not  equal  0  or  1 

Type  IB  -  LUN 

Field  8  not  equal  0,1,2 

Type  1C  -  LUN 

Non-unique  numbers  assigned  to  required  I/O  units. 

Type  1A  -  NJNC 

Field  3  greater  than  44  or  less  than  2. 

Type  1A  -  ANC 

Field  3  greater  than  44  or  less  than  1. 

Type  1A  -  IR 

Field  3  =  Field  4. 

Type  IB  -  IR 

Fields  3  or  4  greater  than  44  or  less  than  1. 

Type  1A  -  CAB 

Field  3  greater  than  22  or  less  than  1. 

Type  IB  -  CAB 

Field  4  =  Field  5. 

Type  1C  -  CAB 

Fields  4  or  5  greater  than  44  or  less  than  1. 

Type  ID  -  CAB 

Fields  7,  8  or  9  less  than  or  equal  to  0. 

Type  IE  -  CAB 

Fields  10  or  11  less  than  0 
Type  IF  -  CAB 

Field  10=0  and  Field  11  not  equal  0 
Type  1G  -  CAB 

Field  10  not  equal  0  and  Field  11=0. 
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Type  1H  -  CAB 

Field  12  greater  than  50  or  less  than  1. 

Type  1A  -  DCAB 

Field  3  greater  than  22  or  less  than  1. 

Type  IB  -  DCAB 

Field  4  greater  than  2  or  less  than  1. 

Type  1C  -  DCAB 

Field  5  greater  than  1000  or  less  than  1. 

Type  ID  -  DCAB 

Field  4=1  and  Field  9  less  than  or  equal  to  0. 

Type  IE  -  DCAB 

Field  4=2  and  Field  9  not  equal  0. 

Type  IF  -  DCAB 

Fields  7,8  or  10  less  than  0. 

Type  1A  -  DJNC 

Field  3  greater  than  44  or  less  than  1. 

Type  IB  -  DJNC 

Field  7  or  8  less  than  0. 

Type  1A  -  DEN 

Field  3  less  than  0. 

Type  1A  -  NDAT 

Field  3  not  equal  0,1,  or  2. 

Type  1A  -  COMP 

Field  3  less  than  or  equal  0. 

Type  1A  -  ANG 

Field  4  less  than  or  equal  0. 

Type  IB  -  ANG 

Field  5  less  than  Field  3. 

Type  2A  -  ANC 

The  junction  number  assigned  to  the  anchor  (Field  3)  has 
been  assigned  to  a  preceding  ANC  card  or  to  an  S=L  junc¬ 
tion  (Field  5)  on  a  preceding  CAB  card. 

Type  2A  -  CAB 

The  junction  number  assigned  to  the  S=L  junction  (Field 
5)  has  been  assigned  on  a  preceding  CAB  card  or  to  an 
anchor  junction  (Field  3)  on  a  preceding  ANC  card. 

Type  3  -  CAB 

A  Type  3  error  appears  only  in  conjunction  with  a  CAB  card 
and  indicates  that  the  number  assigned  to  the  cable  (Field 
3)  has  been  assigned  on  a  preceding  CAB  card. 

Type  4A  -  IR 

The  junction  number  assigned  in  Field  3  has  been  assigned 
in  Fields  3  or  4  of  a  preceding  IR  card. 
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Type  4B  -  IR 

The  junction  number  assigned  in  Field  4  has  been  assigned 
in  Field  3  of  a  preceding  IR  card. 

Type  6A  -  NJNC 

An  NJNC  card  has  previously  appeared  in  the  particular 
source  deck. 

Type  6A  -  DEN 

A  DEN  card  has  previously  appeared  in  the  particular 
source  deck. 

Type  6A  -  COMP 

A  COMP  card  has  previously  appeared  in  the  particular 
source  deck. 

Type  6A  -  VEL 

Twenty-five  VEL  cards  have  previously  appeared  in  the 
particular  source  deck. 

Type  6B  -  VEL 

The  Z-coordinate  at  which  the  current  velocity  is  speci¬ 
fied  (Field  3)  has  been  used  on  a  preceding  VEL  card  in 
the  particular  source  deck. 

Type  7  - 

A  Type  7  error  indicates  an  inadequacy  of  information  in 
cable  array  source  deck  (or  tape) .  The  other  information 
column  under  the  error  heading  contains  a  1x3  matrix,  the 
elements  of  which  give  respectively  the  number  of  NJNC 
cards  read,  the  number  of  DEN  cards  read,  and  the  number 
of  ANC  cards  read.  A  zero  element  is  an  error  (see  cable 
array  source  deck) . 

Type  8  - 

A  Type  8  error  indicates  a  discontinuity  in  numbering 
the  cables  in  the  array.  The  other  information  column 
under  the  error  heading  contains  a  1x22  matrix,  the  ele¬ 
ments  of  which  contain  one  or  zero  indicating,  respec¬ 
tively,  the  use  or  non-use  of  the  corresponding  column 
number  as  a  cable  number.  Zeros  interspersed  with  ones 
are  in  error  (see  array  description). 

Type  9  - 

A  Type  9  error  indicates  a  discontinuity  in  numbering  the 
junctions  in  the  array.  The  other  information  column 
under  the  error  heading  contains  a  1x44  matrix,  the  ele¬ 
ments  of  which  contain  one  or  zero  indicating  respec¬ 
tively,  the  use  or  non-use  of  the  corresponding  column 
number  as  a  junction  number.  Zeros  interspersed  with 
ones  are  in  error  (see  array  description  and  reduction 
to  a  statically  determinate  array). 

Type  11- 

A  Type  11  error  indicates  an  improper  reduction  of  the 
original  cable  array  to  a  statically  determinate  array 
or  an  absence  of  certain  input  cards  in  the  cable  array 
source  deck  (or  tape).  The  other  information  column 
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under  the  error  heading  contains  a  1x5  matrix,  the  ele¬ 
ments  of  which  give,  respectively,  the  number  of  CAB 
cards  read,  the  number  oi  ANC  cards  read,  the  number  of 
junctions  in  the  original  cable  array  (Field  3  of  the 
NJNC  card),  the  number  of  required  cuts  (Eq.  (1))  calcu¬ 
lated  from  the  preceding  information,  and  the  number  of 
IR  cards  read.  Column  5  not  equal  to  column  4  is  an 
error  (see  reduction  to  a  statically  determinate  array). 

Type  12A  -  IR 

The  junction  number  assigned  in  Field  3  is  less  than  or 
equal  to  the  number  of  junctions  in  the  original  (unre¬ 
duced)  array  (Field  3  of  the  NJNC  card). 

Type  12B  -  IR 

The  junction  number  assigned  in  Field  3  is  greater  than 
the  number  of  junctions  in  the  original  (unreduced)  array 
plus  the  number  of  cuts  made  in  reducing  the  array  (Eq. 

(l)). 

Type  12C  -  IR 

The  junction  number  assigned  in  Field  4  is  greater  than 
the  number  of  junctions  in  the  original  (unreduced)  array. 

Type  12A  -  CAB 

The  junction  number  assigned  in  Field  4  is  greater  than 
the  number  of  junctions  in  the  original  (unreduced)  array 
(Field  3  of  the  NJNC  card)  plus  the  number  of  cuts  made 
in  reducing  the  array  (Eq.  (1)). 

Type  12A  -  DCAB 

The  cable  number  assigned  in  Field  3  does  not  correspond 
to  a  cable  number  assigned  to  a  cable  (Field  3  of  the 
CAB  cards). 

Type  12B  -  DCAB 

The  distance  of  the  discrete  device  from  the  S=0  junction 
of  the  cable  (Field  10)  is  greater  than  or  equal  to  the 
length  of  the  corresponding  cable  (Field  9  of  the  CAB 
card) . 

Type  12A  -  DJNC 

The  junction  number  assigned  in  Field  3  does  not  corre¬ 
spond  to  either  a  junction  number  assigned  to  an  anchor 
(Field  3  or  the  ANC  cards)  or  a  junction  number  assigned 
to  an  S=L  junction  of  a  cable  (Field  5  of  the  CAB  cards). 

Type  13  - 

A  Type  13  error  indicates  that  the  original  cable  array 
has  not  been  properly  reduced  to  a  statically  determinate 
array  or  that  junctions  have  been  improperly  numbered. 

The  other  information  column  under  the  error  heading  con¬ 
tains  the  message,  improper  array  reduction  or  junction 
numbering.  Check  tree  representation  of  array  (see 
array  reduction  section  of  Users  Manual)  against  junction 
numbering  on  ANC  and  CAB  cards. 
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Type  14A  - 

A  Type  14  error  indicates  an  inadequacy  of  information 
in  a  parametric  study  source  deck.  The  other  informa¬ 
tion  column  under  the  error  heading  contains  a  1x5  matrix, 
the  elements  of  which  give,  respectively,  the  number  of 
COMP  cards  read,  the  current  field  option,  the  number 
of  VEL  cards  read,  the  number  of  VEL  cards  containing  a 
Z-coordinate  (Field  3)  less  than  or  equal  to  the  minimum 
Z-coordinate  transmitted  by  the  ANC  cards  (Field  6  of 
the  ANC  cards),  and  the  number  of  ANG  cards  read.  If 
column  1  contains  a  zero,  then  a  COMP  card  has  not 
appeared  after  the  first  NDAT  card  (see  parametric 
study  source  decks). 

Type  14B  - 

A  Type  14  error  indicates  an  inadequacy  of  information 
in  a  parametric  study  source  deck.  The  other  information 
column  under  the  error  heading  contains  a  1x5  matrix,  the 
elements  of  which  give,  respectively,  the  number  of  COMP 
cards  less  than  or  equal  to  the  minimum  Z-coordinate 
transmitted  by  the  ANC  cards  (Field  6  of  the  ANC  cards), 
and  the  number  of  ANG  cards  read.  If  column  2  contains 
a  one  and  any  of  columns  3,  4  or  5  contain  a  zero,  then 
the  standard  current  field  has  not  been  properly  formu¬ 
lated  (see  parametric  study  source  decks  and  standard 
current  field). 

Type  15  - 

A  Type  15  error  indicates  than  an  unpermitted  change  has 
been  attempted  in  a  parametric  study  source  deck.  (See 
parametric  study  source  decks.) 

Type  16  - 

A  Type  16  error  indicates  an  improper  deck  structure. 

See  cable  array  source  deck,  parametric  study  source 
decks,  and  overall  input  deck. 

Type  17  - 

A  Type  17  error  indicates  that  the  cable  array  source 
deck  (or  tape)  contains  more  than  2150  records.  The 
other  information  column  under  the  error  heading  contains 
the  message  common/Bl/bounds  exceeded.  See  Users  Manual. 

A  Type  17  error  is  readily  correctable  if  the  machine 
being  used  has  sufficient  core  storage.  This  correction 
is  achieved  by  changing  the  row  dimension  of  DATAT  on 
cards  DES025  and  INF022  from  2150  to  a  number  exceeding 
the  number  of  records  in  the  cable  array  source  deck  (or 
tape).  Simultaneously  the  comparison  value  on  card  INP615 
must  be  changed  from  2150  to  the  new  row  dimension  of 
DATAT. 

Type  18  - 

A  Type  18  error  indicates  that  the  accuracy  required  for 
the  array  equilibrium  calculations  (Field  3  of  the  COMP 
card)  has  not  been  obtained  during  the  calculations. 

This  is  for  one  of  two  possible  reasons. 
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A.  Some  cable  segments  have  gone  slack  (that  is,  the 
segments  have  near  zero  tension.  An  examination  of 
the  tensions  printed  out  in  conjunction  with  a  Type 
18  error  will  reveal  if  this  is  the  reason.  If  it 
is,  see  the  section  on  statically  unstable  cable 
arrays  in  Reference  2  for  possible  remedial  actions. 

If  it  is  not  then  . . . 

B.  The  accuracy  required  for  the  equilibrium  calculations 
is  simply  too  stringent  for  the  computer  to  handle 
(see  COMP  card).  An  examination  of  the  final  value 

of  the  accuracy  obtained,  printed  out  in  conjunction 
with  a  Type  18  error,  will  reveal  the  best  accuracy 
obtainable.  Field  3  of  the  COMP  card  should  be  modi¬ 
fied  to  reflect  this  information. 


ARRAY  DESCRIPTIVE  OUTPUT 

Following  each  error-free  reading  of  a  parametric  study  source 
deck,  DECEL1  transmits  to  the  line  printer  a  description  of  the  physical 
characteristics  of  the  cable  array  under  study.  This  printout  includes 
anchor  junctions  and  locations,  information  concerning  the  reduction  of 
the  original  cable  array  to  a  statically  determine  array,  properties  of 
the  cables  in  the  array,  properties  of  the  discrete  devices  located  at 
array  junctions  total  number  of  Type  1  and  2  devices,  current  field 
option  and  current  profile  and  calculational  accuracy  requirements. 

The  format  of  the  array  descriptive  output  is  self-explanatory.  A 
sample  printout  is  given  in  Appendix  A. 


STRUCTURAL  OUTPUT 

If  output  option  0  or  2  is  selected,  DECEL1  transmits  to  the  line 
printer  a  structural  output.  The  structural  printout  follows  and  refers 
to  the  array  characterized  in  the  array  descriptive  output  and  contains 
information  giving: 

a.  A  description  of  the  current  field  (i.e.,  no  current  or 
current  from  xxx  degrees) 

b.  The  cable  forces  and  angles  at  each  anchor 

c.  The  position  coordinates  of  the  original  (unreduced)  array 
junctions,  the  displacement  of  these  coordinates  from  the 
no-current  coordinates,  and  the  cable  forces  and  angles  at 
each  junction 

d.  The  maximum  and  minimum  tensions  and  their  locations  for 
each  array  cable  and  the  maximum  displacement  from  the  no¬ 
current  condition  and  its  location  for  each  cable 

e.  The  position  coordinates,  the  displacement  of  these  coordi¬ 
nates  from  the  no-current  coordinates,  and  the  tension  at 
each  Type  1  and  2  DCAB  device  in  the  array 


This  latter  information  is  printed  out  in  the  same  order  that  the 
Type  1  and  2  DCAB  devices  are  numbered.  Note  that  it  is  possible  to 
obtain  the  latter  information  for  any  point  on  any  cable  in  the  array 
by  defining  a  "dummy"  Type  2  DCAB  device  to  be  located  at  the  point.  A 
dummy  Type  2  DCAB  device  is  simply  one  for  which  fields  6,  7,  and  8  of 
the  DCAB  card  are  left  blank  so  that  the  "device"  has  no  effect  on  the 
array  equilibrium  calculations. 

The  format  of  the  structural  output  is  self-explanatory.  A  sample 
printout  is  given  in  Appendix  A  (Figure  17). 


DEVICE  LOCATION  OUTPUT 

If  output  option  1  or  2  is  selected,  DECEL1  transmits  to  tape  or 
cards  a  device  location  output.  The  device  location  output  contains 
information  giving: 

1.  A  description  of  the  current  field  (i.e.,  no  current  or 
current  from  xxx  degrees) 

2.  The  position  coordinates  of  each  Type  1  and  2  discrete 
device  in  the  array 

Four  types  of  records  are  associated  with  the  device  location 
output.  Each  of  these  records  is  written  with  the  format: 

(A4,I4,3F10.2) 
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REC  Record 


Field 

Format 

Contents 

1 

A4 

.REC 

2 

14 

Record  number 

3 

F10.2 

Not  used 

4 

F10.2 

Not  used 

5 

F10.2 

Not  used 

NOTE:  The  REC  record  is  used  to  identify  the  cable 

array  under  study.  The  record  number  (Field  2) 
is  referenced  in  the  array  descriptive  output 
in  the  statement,  "DEVICE  LOCATION  OUTPUT 
RECORD  XX  REFERS  TO  THIS  ARRAY."  The  informa¬ 
tion  following  a  REC  record  refers  to  the  iden¬ 
tified  array. 


CUR  Record 


Field 

Format 

Contents 

1 

A4 

.CUR 

2 

14 

0  if  no  current  acting  on  array 

1  if  current  is  acting  on  array 

3 

F10.2 

Blank  if  Field  2=0 

Current  angle  if  Field  2=1 

4 

F10.2 

Not  used 

5 

F10.2 

Not  used 

NOTE:  The  CUR  record  is  used  to  describe  the  current 
field  (i.e.,  no  current  or  current  from  xxx 
degrees).  The  information  following  a  CUR 
record  refers  to  the  identified  current  field. 
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DEV  Record 


Fi  -'Id 

Format 

Contents 

1 

A4 

.DEV 

2 

14 

Discrete  device  index  (Field  5 
of  the  DCAB  and  DJNC  cards) 

3 

F10.2 

X  coordinate  of  device 

4 

F10.2 

Y  coordinate  of  device 

5 

F10.2 

Z  coordinate  of  device 

NOTE:  The  DEV  records  are  used  to  transmit  the  posi¬ 
tion  coordinates  of  the  Type  1  and  2  discrete 
devices  in  the  array.  There  is  one  DEV  record 
for  each  indexed  (Types  1  and  2)  device  for 
each  identified  current  field. 


EOP  Record 


Field 

Format 

Contents 

1 

A4 

.EOP 

2 

14 

Not  used 

3 

F10.2 

Not  used 

4 

F10.2 

Not  used 

5 

F10.2 

Not  used 

NOTE:  The  EOP  record  is  used  to  signify  the  end  of 
output  transmission. 
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A  typical  overall  device  location  output  file  is  illustrated  below. 

-REC-^1 

~CUR„„..0 

DEV  Records  (one  for  each  indexed  device) 

-CUR.,,,.1  current  angle 
DEV  Records 

..CUR.. ....  1  current  angle 

DEV  Records 
..REC......2 

-.CUR......0 

DEV  Records 

..EOP 
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Appendix  A 
EXAMPLE  PROBLEM 
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An  example  problem  has  been  included  both  to  illustrate  actual 
input  and  output  and  as  a  test  case.  The  test  case  can  be  used  to 
confirm  proper  operation  of  DECEL1  when  implemented  on  a  particular  host 
computer.  The  test  case  structure  is  shown  in  Figure  13.  The  structure 
represents  an  acoustic  array  in  the  horizontal  leg;  this  leg  has  been 
buoyed  at  the  center  to  keep  it  approximately  horizontal.  A  signal 
cable  rises  from  one  anchor  to  a  subsurface  buoy.  The  applied  current 
is  unidirectional  with  the  profile  shown  in  Figure  14.  The  current  is 
acting  broadside  to  the  structure.  The  DECEL1  model  of  the  structure  is 
shown  in  Figure  15;  other  details  of  the  modeling  are  shown  in  Tables  4 
through  8. 

Input  card  images  are  shown  in  Figure  16.  Note  that  one  dummy-DCAB- 
device  card  was  intentionally  placed  "out  of  sequence"  just  after  the 
NJNC  card.  This  illustrates  that  the  card  sequence  is  essentially 
arbitrary,  except  as  noted  elsewhere  in  the  manual. 

The  output  from  DECEL1  is  shown  in  Figures  17  and  18.  Note  that 
the  "number  of  indexed  devices"  count  includes  the  three  DJNC  devices. 

The  "Array  Equilibrium  with  no  Current"  portion  defines  both  the  initial 
positions  of  devices  and  the  properties  of  the  DCAB  devices.  If  the  "do 
not  print"  flag  had  been  set  on  any  DCAB  card,  that  device  would  not 
have  been  tabulated  here;  however,  its  effects  would  have  been  accounted 
for  in  the  solution. 

In  the  portion  of  the  output  where  the  current  was  applied  (Figure 
18),  the  displacement  of  all  junctions  and  devices  is  listed  relative  to 
the  present  no-current  position  and  relative  to  the  original  no-current 
case.  In  this  case  both  displacements  are  identical  since  the  new 
no-current  reference  flag  had  not  been  set  on  the  NDAT  card. 

Figure  19  shows  the  configuration  of  the  array  both  in  the  no-current 
condition  (dotted  lines)  and  with  current  applied  (solid  lines).  For 
this  plot  the  default  perspective  view  angles  were  used.  Other  views 
(plan  or  elevation)  of  the  same  structure  in  the  same  current  could  be 
obtained  by  additional  identical  NDAT  cases  with  the  perspective  plotting 
angles  changed  accordingly. 


80 


Table  4.  Cables 


Cable 

No. 

Junction 

Length 

(ft) 

Weight/ 

Foot 

(lb/ft) 

Diameter 
(in. ) 

Drag 

Coefficient 

From 

To 

1 

1 

2 

3,000 

0.25 

1.0 

1.4 

2 

2 

8 

3,000 

0.25 

1.0 

1.4 

3 

2 

4 

5,000 

0.20 

1.4 

4 

4 

9 

3,000 

0.25 

1.0 

1.4 

5 

4 

10 

3,000 

0.25 

1.0 

1.4 

6 

9 

7 

4,000 

0.3 

1.25 

1.4 

Table  5.  Anchor  Locations 


Anchor  No. 

Junction  No. 

X 

y 

z 

1 

1 

0 

1,000 

0 

2 

3 

0 

-1,000 

0 

3 

5 

7,000 

-1,000 

0 

4 

6 

_ 

7,000 

1,000 

0 

Table  6.  Imaginary  Reaction  Cuts 


Junction  No. 

Cut  No. 

3 

8 

5 

9 

6 

10 

Note:  In  the  model  the  cables  terminate  at 
cuts  rather  than  at  anchors. 
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Table  7.  DJNC  Devices 


Device 

Junction 

No. 

Device 

Buoyancy 

(lb) 

Device  Drag 
Coefficient 

Device 

Frontal 

Area 

(ft2) 

2 

3,000 

1.0 

26.0 

4 

3,000 

1.0 

26.0 

7 

3,000 

1.0 

26.0 

Table  8.  DCAB  Devices 


Device  Index 

On  Cable  No. 

S  Coordinate 
(ft) 

Device 

Buoyancy 

(lb) 

9 

3 

2,500 

1,000 

Note:  Other  dummy  devices  are  used  to  obtain  a  printout  of  the 
spatial  location  of  the  device;  these  devices  have  not 
been  tabulated  here.  Dummy  devices  have  a  very  small 
weight  and  no  drag. 
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Figure  16.  Input  card  images. 
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Figure  17.  Continued. 
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AFB  CESCH,  Wright-Patterson 

ARMY  ARRADCOM.  Dover.  NJ 

ARMY  -  CERL  Library,  Champaign  IL 

ARMY  COASTAL  ENGR  RSCH  CEN  Fort  Belvoir  VA 

ARMY  COE  Philadelphia  Dist  (LIBRARY)  Philadelphia.  PA 

ARMY  ENG  WATERWAYS  EXP  STA  Library.  Vicksburg  MS 

ARMY  ENGR  DIST.  Library.  Portland  OR 

ARMY  MOBIL  EQUIP  R&D  COM  Fuel  Hndlng  Equip  Br..  Ft  Belvoir.  VA;  Mr.  Ccvasco.  Fort  Belvoir  MD 
CNM  MAT  08T245  (Spalding).  Washington.  DC 
COMSUBDEVGRUONE  Operations  Offr.  San  Diego,  CA 
DOE  Dr.  Cohen 

DTNSRDC  Code  522  (Library).  Annapolis  MD;  Pattison.  Code  1706.  Bethesda,  MD;  Rispin,  Code  1706. 
Bethesda,  MD 

LIBRARY  OF  CONGRESS  Washington.  DC  (Sciences  &  Tech  Div) 

NATL  RESEARCH  COUNCIL  Naval  Studies  Board.  Washington  DC 
NAVCOASTSYSTCTR  Library  Panama  City.  FL 

NAVELEXSYSCOM  Code  PME-124-61.  Washington  DC;  PME  124-30  (Henderson)  Washington.  DC 
NAVFACENGCOM  Code  0453  (D.  Potter)  Alexandria.  VA;  PC-2  Alexandria,  VA 
NAVFACENGCOM  -  CHES  DIV.  Code  FPO-1E.  Wash.  DC 
NAVOCEANO  Code  1600  Bay  St.  Louis.  MS 

NAVOCEANSYSCEN  Code  4473  Bayside  Library,  San  Diego.  CA;  Code  5204  (J.  Stachiw),  San  Diego,  CA; 

Tech.  Library.  Code  447 
NAVSCOLCECOFF  C35  Port  Hueneme.  CA 
NAVSEASYSCOM  Code  SEA  OOC  Washington,  DC 
NAVSEC  Code  6034  (Library).  Washington  DC 
NAVSHIPREPFAC  Library.  Guam 

NAVSHIPYD  Code  202.4,  Long  Beach  CA;  Code  202.5  (Library)  Puget  Sound.  Bremerton  WA;  Library, 
Portsmouth  NH;  Tech  Library.  Vallejo.  CA 
NOAA  Library  Rockville,  MD 

NOAA  DATA  BUOY  OFFICE  Engrng  Div  (Riannie)  Bay  St.  1-ouis,  MS 
NORDA  CO,  Bay  St.  Louis,  MS;  Code  350  (Swenson)  Bay  St.  Louis,  MS 

NRL  Code  2627.  Washington  DC;  Code  8400  Washington.  DC;  Code  8441  (R.A.  Skop).  Washington  DC 
NUSC  Code  131  New  London,  CT;  Code  S332,  B-80  (J.  Wilcox) 

OCEANAV  Mangmt  Info  Div..  Arlington  VA 

ONR  (Dr.  E.A.  Silva)  Arlington.  VA;  Central  Regional  Office.  Boston.  MA 

PMTC  Code  4253-3,  Point  Mugu,  CA 

PWC  Code  120.  Oakland  CA;  Code  154.  Great  Lakes.  IL 

UCT  TWO  OIC,  Norfolk.  VA 

U  S.  MERCHANT  MARINE  ACADEMY  Kings  Point.  NY  (Reprint  Custodian) 

US  NATIONAL  MARINE  FISHERIES  SERVICE  Highlands  NY  (Sandy  Hook  Lab-Library) 

USCG  (G-MP-3/USP/82)  Washington  Dc 

USCG  R&D  CENTER  CO  Groton,  CT;  D.  Motherway,  Groton  CT;  D.  Paskausky.  Groton.  CT;  Tech.  Dir. 
Groton.  CT 

USNA  Ch.  Mech.  Engr.  Dept  Annapolis  MD;  Ocean  Sys.  Eng  Dept  (Dr.  Monney)  Annapolis.  MD 
AMERICAN  CONCRETE  INSTITUTE  Detroit  Ml  (Library) 

CALIF.  MARITIME  ACADEMY  Vallejo.  CA  (Library) 

CALIFORNIA  INSTITUTE  OF  TECHNOLOGY  Pasadena  CA  (Keck  Ref  Rm) 

CATHOLIC  UNIV.  Mech  Engr  Dept.  Prof.  Niedzwecki.  Wash.,  DC 
CORNELL  UNIVERSITY  Ithaca  NY  (Serials  Dept.  Engr  Lib.) 

DAMES  &  MOORE  LIBRARY  LOS  ANGELES.  CA 
DUKE  UNIV  MEDICAL  CENTER  B  Muga,  Durham  NC 
FLORIDA  ATLANTIC  UNIVERSITY  Boca  Raton  FL  (W.  Tessin) 

LEHIGH  UNIVERSITY  Bethlehem  PA  (Linderman  Lib.  No.3<).  Flecksleiner) 

MAINE  MARITIME  ACADEMY  CAST1NE.  ME  (LIBRARY) 
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MIT  Cambridge  MA;  Cambridge  MA  (Rm  10-500,  Tech.  Reports,  Engr.  Lib.) 

NATL  ACADEMY  OF  ENG.  ALEXANDRIA,  VA  (SEARLE,  JR.) 

OREGON  STATE  UNIVERSITY  (CE  Dept  Grace)  Corvallis,  OR 
PENNSYLVANIA  STATE  UNIVERSITY  State  College  PA  (Applied  Rsch  Lab) 

PURDUE  UNIVERSITY  Lafayette,  IN  (CE  Engr.  Lib) 

SOUTHWEST  RSCH  INST  R.  DeHart,  San  Antonio  TX 
STANFORD  UNIVERSITY  Engr  Lib,  Stanford  CA 
STATE  UNIV.  OF  NEW  YORK  Buffalo,  NY 

UNIVERSITY  OF  HAWAII  HONOLULU,  HI  (SCIENCE  AND  TECH.  DIV.);  Ocean  Engmg  Dept 
UNIVERSITY  OF  ILLINOIS  URBANA,  IL  (LIBRARY) 

UNIVERSITY  OF  MASSACHUSETTS  (Heronemus),  Amherst  MA  CE  Dept 
UNIVERSITY  OF  NEW  HAMPSHIRE  (Corell)  Durham,  NH 
UNIVERSITY  OF  NOTRE  DAME  Katona,  Notre  Dame,  IN 
UNIVERSITY  OF  RHODE  ISLAND  Narragansett  RI  (Pell  Marine  Sri.  Lib.) 

UNIVERSITY  OF  SO.  CALIFORNIA  Univ  So.  Calif 
UNIVERSITY  OF  TEXAS  Inst.  Marine  Sri  (Library),  Port  Arkansas  TX 
AMETEK  Offshore  Res.  &  Engr  Div 

COLUMBIA  GULF  TRANSMISSION  CO.  HOUSTON,  TX  (ENG.  LIB.) 

EG&G  WASH  ANALYTICAL  SERV  CTR,  INC  Rockville,  MD 
EXXON  PRODUCTION  RESEARCH  CO  Houston  TX  (A.  Butler  Jr) 

GOULD  INC.  Shady  Side  MD  (Ches.  Inst.  Div.,  W.  Paul) 

GRUMMAN  AEROSPACE  CORP.  Bethpage  NY  (Tech.  Info.  Ctr) 

HONEYWELL,  INC.  Minneapolis  MN  (Residential  Engr  Lib.) 

MOBIL  PIPE  LINE  CO.  DALLAS,  TX  MGR  OF  ENGR  (NOACK) 

OCEAN  ELECTRONICS  APPLICATION,  INC.  (Softley)  Key  Biscayne,  FL 
SANDIA  LABORATORIES  Library  Div.,  Livermore  CA 
TRW  SYSTEMS  CLEVELAND,  OH  (ENG.  LIB.) 

WESTINGHOUSE  ELECTRIC  CORP.  Annapolis  MD  (Oceanic  Div  Lib,  Bryan) 

WESTINTRUCORP  Egerton,  Oxnard,  CA 

WM  CLAPP  LABS  -  BATTELLE  DUXBURY,  MA  (LIBRARY) 
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